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Abstract — In this paper, we consider CSMA policies for 
scheduling of multihop wireless networks with one-hop traffic. 
The main contribution of this paper is to propose Unlocking 
CSMA (U-CSMA) policy that enables to obtain high throughput 
with low (average) packet delay for large wireless networks. 
In particular, the delay under U-CSMA policy becomes order- 
optimal. For one-hop traffic, delay is defined to be order-optimal 
if it is O(l), i.e., it stays bounded, as the network-size increases to 
infinity. Using mean field theory techniques, we analytically show 
that for torus (grid-like) interference topologies with one-hop 
traffic, to achieve a network load of p, the delay under U-CSMA 
policy becomes 0(1/(1 p) 3 ) as the network-size increases, and 
hence, delay becomes order optimal. We conduct simulations 
for general random geometric interference topologies under U- 
CSMA policy combined with congestion control to maximize 
a network-wide utility. These simulations confirm that order 
optimality holds, and that we can use U-CSMA policy jointly 
with congestion control to operate close to the optimal utility 
with a low packet delay in arbitrarily large random geometric 
topologies. To the best of our knowledge, it is for the first time that 
a simple distributed scheduling policy is proposed that in addition 
to throughput/utility-optimality exhibits delay order-optimality. 

I. Introduction 

One of the most intriguing challenges in the context of 
wireless networking is the design of a scheduling policy that 
has the following properties: 

a) throughput-optimality, 

b) low packet delajQ and 

c) simple and fully distributed implementation. 

From a complexity theoretic viewpoint, unless NP C BPP 
or P = NP, there does not exist [1] a universal scheduling 
policy that has the above three properties for all possible net- 
work topologies. However, it is still possible to design a policy 
that has the above properties for a subset class of network 
topologies. This seems to be true for geometric networks [2], 
[3], in which only links that are geometrically close interfere 
with each other. These networks closely approximate a wide 
range of practical wireless networks, and yet are known to 
admit Polynomial-Time Approximation Scheme (PTAS) for 
several NP-hard optimization problems (see e.g., [2], [4]). For 
this reason, our focus in this paper is to design a scheduling 
policy for large geometric wireless networks. 

There are two main approaches into the design of scheduling 
policies in wireless networking: either through matching poli- 

1 Throughout this paper, by delay we mean average packet delay. 



cies [2], [3], [5]— [16], or through random access policies [17]- 
[36]. Despite the past efforts that have significantly advanced 
our understanding of these policies and their performance, to 
the best of our knowledge, there is no instance of these policies 
that realizes all of the three properties mentioned earlier, even 
for geometric networks. 

On one hand, we have matching policies that can be 
throughput-optimal [5]-[7], [9] and can provide order-optimal 
low delay [15]. However, these optimalities are obtained 
assuming that an NP-hard problem can be solved in each 
scheduling round. At the same time, reducing the complexity 
of matching policies, in general, comes at the price of losing 
throughput-optimality of these policies [10]— [13] or a large 
delay [1]. This leaves the design of a matching policy with 
all the three properties as an open research challenge (see 
Section [It] for further discussion on matching policies). 

On the other hand, we have random access policies that 
are naturally simple and can be implemented in distributed 
manner. Among these policies, the classical CSMA policy, 
see [24], [27], [29]— [3 1] for variants of this policy, is the 
one that features throughput-optimality in a wide range of 
settings [24], [29]-[31]. This encourages to use the classical 
CSMA policy for scheduling of wireless networks. However, 
as will be discussed shortly, the delay performance under this 
policy can be very poor. As a result, the current random access 
policies do not possess all of the three properties mentioned 
earlier. 

As a motivating example, consider annxn torus (see Fig.[T| 
interference graph [2], [29]-[31] with L — n 2 nodes where 
each node interferes with the four closest neighbouring nodes. 
Suppose packet arrival rate is uniform, i.e., it is the same for 
all nodes, and that is equal to A. Let p be the corresponding 
loacQ and define 

e = 1 - p. 

For this simple topology, a mixing-time analysis [31] upper- 
bounds the packet delay under the classical CSMA policy as 



O 



where c„ > 1 is a constant. For small e, a similar analysis 
[37] lowerbounds the packet delay under the classical CSMA 

2 In the limit of large toruses, the maximum unifor m throughput is 0.5, and 



load p in the limit becomes See Section III-C for the definition of p. 



policy as 



,ciL/(logLY 



), 



for some constant c; > 0. 

The above delay-bounds show that the classical CSMA 
policy exhibits a threshold behaviour in the sense in order 
to achieve a high throughput, i.e., to make e small, one has 
to tolerate a delay that exponentially grows with the network- 
size L. The threshold behaviour and the exponential growth are 
related to the phase transition phenomenon^] in the hard-core 
lattice gas model [38], [39]. Due to such threshold behaviours, 
even in mid-sized simple topologies, the classical CSMA 
policy cannot support a high throughput with low delay (see 
Section [rV-A) i. 

In this paper, we propose Unlocking CSMA (U-CSMA) as 
a new CSMA policy that overcomes the threshold behaviour 
of the classical CSMA policy. While being simple and dis- 
tributed, U-CSMA policy has the following properties for 
geometric networks with one-hop traffic [2], [3]. 

a) It enables to achieve a high throughout/utility arbitrarily 
close to the optimal with a low (average) packet delay. 

b) The (average) packet delay under this policy is order- 
optimal, i.e., it stays bounded as the network-size in- 
creases to infinity. 

We provide analytical results for the torus interference 
topology with uniform packet arrival rate as considered earlier, 
and show that for large network-size L, the average delay 
under U-CSMA policy is order-optimal and is 



o 



(1) 



It is important to note that the above delay bound is indepen- 
dent of the network size L, in sheer contrast to the delay under 
the classical CSMA policy that exponentially increases with 
the network-size L. This means that U-CSMA policy does 
not suffer from the threshold behaviour and is indeed able to 
provide high throughput with low delay for arbitrarily large 
torus topologies. 

In our simulation study, we use U-CSMA policy jointly 
with a congestion control algorithm to maximize a network- 
wide utility in large random geometric networks. We show that 
using U-CSMA, we can assign packet arrival rates closely to 
the optimal with a low packet delay that stays bounded as 
the network-size increases, and hence, a delay that exhibits 
order-optimality. As far as we are aware, it is for the first 
time that a simple distributed scheduling policy is proposed 
that can operate close to the optimal with order-optimal low 
packet delay. 

We believe that the design principle of U-CSMA policy and 
the novel approach taken to study its performance open up a 
new direction into the design and study of scheduling policies 
for large-scale wireless networks. The main significance of our 
study in this paper is that it realizes the possibility of having 
large-scale wireless multihop networks that can be maintained 

3 Phase transition has also been reported as the cause of border effects that 
persist in 2D under the classical CSMA policy [27]. 



in a simple distributed manner and that can provide high 
throughput/utility, arbitrarily close to the optimal, with order- 
optimal low packet delay. 

A key step to obtain the delay bound in ([T| is where 
we show that the schedule under the classical CSMA pol- 
icy quickly converges to a maximum schedule in geometric 
networks. Using techniques from mean field theory [40], we 
show that for large torus and lattice topologies with large 



uniform attempt-rates, the distance (see Section VI-Ai to the 



maximum schedules as a function of time t drops as . To the 
best of our knowledge, our result is the first that analytically 
characterizes the fast convergence behaviour of the classical 
CSMA policy. As this convergence is independent of network- 
size L, it is fundamentally different than the convergence time 
to the steady-state (i.e., the mixing time) of the dynamics of 
the classical CSMA policy, which can be exponentially large 
in L [37]. 

The rest of the paper is organized as follows. In the next 
section, we briefly review the related work. In Section [Hi] 
we present the network model and the classical CSMA policy 



model. In Section IV we provide an overview of our main 
results, including the description of U-CSMA policy and sim- 
ulation results. In Section[V] we provide one example to imple- 
ment U-CSMA policy in a fully distributed and asynchronous 
manner. In section [Vl] we provide a formal statement of our 



analytical results in this paper. In Section VII we elaborate on 



the dynamics of schedules under the classical CSMA policy 
whose characterization is required to derive analytical results. 
In Section [VIII| we formally state two assumptions that allow 
the formal analysis developed in this paper. In Section [rX] we 
provide details on how to use U-CSMA policy jointly with a 
congestion control algorithm for general topologies. Finally, 
we conclude the paper in Section [X] 

II. Related Work 

In this section, we provide a brief, by no means exhaustive, 
overview of the work in the area of wireless scheduling that 
is closest to ours in this paper. We consider two main classes, 
i.e., the matching policies and random access policies. 

Matching Policies: Maximum Weight Matching (MWM) 
policy was first proposed in the seminal work in [5]. This 
policy is perhaps the first policy that is throughput-optimal 
in a wide range of settings [5]-[7], [9]. MWM policy at any 
timeslot maximizes a weighted summation of queue-sizes in 
the network, which can be an NP-hard optimization problem 
[2]. Despite its complexity, simulations [16] show that MWM 
policy is close to the optimal in terms of delay for one-hop 
traffic. For multihop traffic, the delay under MWM policy is 
0(f), and for one-hop traffic is order-optimal as 0(-), under 
certain conditions [15] that hold for geometric networks. The 
delay bound in our paper for one-hop traffic is 0([^] 3 ), which 
includes a multiplicative factor of [-] 2 as well as -. This 
factor can be interpreted as the scheduling-time needed to find 
schedules that are e close to the optimality. However, we note 
that the delay performance in [16] and the O(-) bound in [15] 



are obtained assuming that the NP-hard problem of MWM 
policy can be solved at every timeslot. 

Greedy Maximal Matching (GMM) policy is a simple and 
distributed alternative for MWM policy, see e.g., [2], [11]. 
While GMM policy is not throughput-optimal in general, a 
number of local pooling results [3], [10], [13] indicate that 
for a noticeable subset of topologies, GMM policy is indeed 
throughput-optimal. However, GMM requires message pass- 
ing, and it is an open area to investigate the delay performance 
of GMM policy. Maximal Matching (MM) policy is simpler 
than GMM policy and has order-optimal delay of O(^) for 
one-hop traffic [12]. However, this policy is not throughput- 
optimal and is guaranteed to stabilize only half of the capacity 
region. In [41], a matching policy is proposed that can stabilize 
arbitrarily close to 100% of the capacity region in expense 
of increasing an overhead that is constant in network-size. 
However, this policy is limited for networks with primary 
interference. See [14] for a comparison of different matching 
policies. 

Random Access Policies: Random access policies started 
with the classical Aloha protocol [17], for which an optimality 
result was first established in [18]. The capacity of random 
access policies under collision detections, acknowledgements, 
or backoff schemes have been studied in [19], [20], [22]. The 
recent work in [26] chooses access probabilities in an Aloha- 
like policy based on queue backlogs to achieve the capacity 
region of slotted Aloha. In [25], [33], distributed protocols 
are proposed that assign access probabilities to maximize a 
network utility under an Aloha-like protocol. Due to their 
simplicity, Aloha-like protocols have been also used in mobile 
networks [23]. These protocols however are not throughput- 
optimal [26]. 

CSMA policies are a special class of random access poli- 
cies that assume nodes can sense whether their neighbours 
are transmitting. Performance of these policies as defined 
in 802.11 standard for a specific network setup is studied 
in [21]. For an interesting but special class of networks 
with primary interference, it is known that 1) CSMA polices 
are throughput-optimal [24], and 2) for a subclass of these 
networks such as the n x n switch, the delay to access the 
channel becomes memoryless under CSMA policies, leading 
to an 0(\) (normalized) packet delay [32]. 

Throughput-optimality of CSMA policies extends to net- 
works with arbitrary interference graphs [29]-[31]. The 
throughput-optimal CSMA policies in [29]-[31] are based on 
a continuous time Markov chain that prevents collisions. This 
is addressed by considering contention resolution [30], [34]. 

Both in [29] and [30], it is assumed that there is a time-scale 
separation and, hence, CSMA dynamics quickly converges to 
its steady-state faster than the rate by which queues change 
over time. The authors of [31] and later those of [36] show that 
as long as attempt rates of nodes change sufficiently slowly, 
throughput optimality can be achieved. A related work [35] 
divides the time axis into frames, and updates parameters of 
CSMA policy only at the beginning of each frame. How- 
ever, delay performance under the above throughput-optimal 



schemes is not investigated, and the upperbound on the delay 
inferred from these papers increases with the network-size. 

Before concluding this section, we note that there are nu- 
merous results that study link starvation under CSMA policies, 
e.g., see [28] and references therein. In particular, the work in 
[27] shows that in 2D, the phase transition phenomenon makes 
the CSMA policy lock into a certain similar set of states for a 
long time, causing large packet delays. Using this insight, we 
propose U-CSMA policy that benefits from a novel unlocking 
mechanism. In cotrast to previous matching or random access 
policies, U-CSMA is a simple and distributed policy that 
provides high throughput with low delay that features order- 
optimality. 

III. Network and Classical CSMA Policy Model 

In this section, we introduce the network and classical 
CSMA policy model that we use in this paper. 

A. Network Model 

We consider a fixed wireless network consisted of a set 
Af of nodes, and a set £ of links with cardinality L. We 
refer to L as the network size. A link I = (n, m) € £ 
indicates that transmitter node n and receiver node m are 
within transmission range of each other and can exchange data 
packets. Each link / = (n, to) corresponds to a queue that is 
maintained by its transmitter node n. 

We model the contention between links by an interference 
graph G{£, £) [2], [29]-[31], [35], where £ is the set of links 
and £ is the set of edges. An edge e = (I, I') e £ in the graph 
G(£,£) indicates that the two links I and I', 1,1' £ £ interfere 
with each other. In the following, we will refer to £ as the 
node set of the interference graph, and to the set £ as its edge 
set. We define a geometric interference graph [2]-[4] to be a 
graph whose vertices can be considered as points on the plane, 
and where two vertices are connected by an edge if and only if 
the distance between them is less than the interference range r 
where r > 0. We define a geometric network as a network with 
geometric interference graph. We define a random geometric 
network as a geometric network for which the vertices of its 
interference graph are points that are distributed according to 
a uniform stochastic process over a convex region in the plane. 

We define a valid schedule to be a subset of links in £ no 
two of which interfere with each other. We define a maximum 
schedule to be a valid schedule with the largest number of 
links in £. We also define a link to be active at time t, if the 
link is transmitting at time t. We define a scheduling policy to 
be an algorithm, randomized or deterministic, that determines 
which links are active at any given time. 

Throughout the paper, we assume that traffic is one-hop. 
Let A; be the packet arrival rate for transmission over link I, 
which corresponds to a queue in the network, and let 

be the arrival rate vector for a given network. We assume that 
the rate of transmission is the same for all links, and it takes 
one unit of time to transmit any one packet. 



To characterize the arrival process in further detail, for t 2 > 
t\ > 0, let Ai(ti,t2) be the number of packets that arrive for 
transmission to link I in the time interval (ti,t?\. We assume 
that the number of packets that arrive in a unit time interval 
to any link I is bounded by a constant A max , i.e., 



Mt,t+1) < A v 



(2) 



Moreover, we assume for any e > 0, there exists an integer 
k e > 1 such that for t% — 1\ > k e and for all / E C, we have 



E 



Ai{t x ,h) 
t 2 -ti 



U s (h) 



< e 



(3) 



where H s (t±) is the system history up to and including time 
t\. The above intuitively means that the expected time-average 
number of packets that arrive to a link I converges to its arrival 
rate A;. 

B. Classical CSMA Policy 

For our analysis, we define the classical CSMA policy as 
follows, similar to the ones presented in [27], [29]— [31]. Given 
a wireless network with interference graph G (£,£), every 
link I £ C independently of others senses transmissions of 
any conflicting link in the interference graph G(C,£), i.e. 
of any link I' such that the edge e = (I, I') is contained in 
the edge set £. A link I senses the channel as idle at time t 
if all of its conflicting (interfering) links are not active and 
not transmitting at time t. If link I senses that any of its 
interfering links is transmitting, then it waits until all of its 
interfering links become silent. Once this happens, link I sets 
a backoff timer with a value that is exponentially distributed 
with mean Zi > 0, and starts to reduce the backoff 

timer. If the timer reaches zero before any of its interfering 
links start a transmission, then link I starts a transmission. 
Otherwise, link I simply waits until all of its interfering links 
become silent again, and repeats the above process. We define 
Zi to be the transmission attempt-rate of link I. We assume 
that all transmission times are independently and exponentially 
distributed with unit mean. 

The above models an idealized CSMA policy in which 
1) any link can always sense transmissions of all of its 
interfering links, and 2) there is no hidden-terminal problem 
that can create packet collisions as in [27], [29], [31]. These 
assumptions can be removed using the methods of [30], [34]. 
Hence, we continue assuming that the above two assumptions 
hold. 

We characterize a classical CSMA policy by the vector z = 
{zi)i£C where z; is the transmission attempt-rate of link I. 
Given vector z, the network dynamics as which links are active 
over time can be represented by a Markov process [29]. Using 
this, we can define /i;(z), I e C, as the service rate of link I 
under z, i.e., /x;(z) is the fraction of time that link / is active 
under the CSMA policy z. 

We say that the classical CSMA policy z stabilizes the 
network for a given packet arrival rate vector A if [5] 



(n,n) 
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A; < m(z), l e C. 



(4) 



Fig. 1. Lattice and torus interference graphs. Each dark circle represents a 
link, and an edge between two dark circles shows that their corresponding 
links interfere with each other. 

This commonly used stability criteria [5] requires that for each 
link I £ C, the link service rate /i; (z) is larger than the arrival 
rate A;. Given a fixed network, we then define the achievable 
rate region C of the classical CSMA policy as 

C = {A : 3z s.t. Q holds.}, 

i.e., as the set of all rate vectors A for which there exists a 
vector z that stabilizes the network for A. 

It is well-known that the classical CSMA policy is through- 
put optimal [24], [29]— [31], i.e., the set C contains all arrival 
rate vectors A that are inside the capacity region T, where T 
is the set of all A's that can be stabilized by any scheduling 
policy, CSMA or not, including those with the full knowledge 
of future packet arrivals. 

C. Lattice and Torus Interference Graphs with Uniform At- 
tempt and Packet Arrival Rates 

To obtain analytical results, we consider wireless networks 
with grid-like interference graphs. In particular, we consider 
the lattice interference graph Gl = Gl{£,£) and the torus 
interference graph Tl = 71(£, £7-)- I n both cases, the set £ is 
the set of all links where each link I £ C can be represented by 
coordinates i,j G {0, n}, on the plane. See Fig.[TJfor 

an illustration. Hence, the network-size, i.e., the total number 
of links, is given by L — (n + l) 2 . 

It remains to specify which links interfere with each other. 
For the lattice interference graph Gl, we assume that there 
exists an edge e 6 £ between any two links I = and 
V = (i',f), 1,1' G C, iff link I and link I' differ in exactly 
one coordinate, i.e., we have that 

|i-i'| + |i-j'| = l. 

For the torus interference graph 71, the edge set £j- contains 
all edges defined for the lattice interference graph G^. In 
addition, the set £7- contains an edge between link I — (i, 0) 
and link /' = (i,n), for < i < n, and also contains an edge 
between link / = (0,j) and link /' = (n,j), for < j < n. 
As a result, the torus interference graph 71 is the same as 
Gl with additional edges around the boundary of Gl so that 
every link has exactly four interfering links. 

Given a lattice or torus interference graph, we define a link 
I = G C as an even link iff i + j is an even number. 

We define £( e ' as the set of all such even links. Similarly, we 



define a link I — (i, j) £ C as an odd link iff i + j is an odd 
number, and define £W as the set of all odd links. 

For the lattice and torus interference graphs Gl and Tl, 
we focus on CSMA policies {z} with uniform transmission 
attempt-rates so that 



zi 



I e C, 



for some z > 0. In addition, we focus on the case of uniform 
packet arrival rates, i.e., we let 



A; = A, 



< A < 

f-^rnax 

(L), I G C. 



(5) 



where p m ax(L) is the maximum uniform-throughput, i.e., the 
maximum throughput that can be provided for all links by any 
policy in the network. For lattice interference graph Gl, we 
have that 

pmax{L) = 0.5. 

This throughput can be achieved, for instance, by alternating 
between two valid schedules £(°) and every unit of time, 
which allows every link to be active half of the time. For 
torus interference graph Tl, due to boundaries being wrapped 
around, £(°) and are not valid schedules, but we can show 
that 

lim p max {L) = 0.5. 

Having defined p max (L), for a given lattice or torus inter- 
ference graph with L links, we define the network load factor 
or simply load p as 

p = pW = — Kr v (6) 

f-^max \-'- J ) 

We also define e to be the distance to maximum load of p = 1: 
e = e (A) = 1 - p(A). (7) 
We next provide an overview of our main results. 

IV. Overview of Main Results 

In this section, we provide an overview of our main results. 
We first investigate the performance of the classical CSMA 



policy as defined in Section III-B and explain why under this 
policy it is impractical to obtain both high throughput and 
low delay. We explain that a locking-in behaviour leads to an 
exponentially increasing delay that for the torus interference 
graph is upperbounded as 



O 



(IT) 



where e is defined in Q, and c u > 1 is a constant. We also 
explain that for small e, the delay is 



ciL/ (log L) 



), 



for some constant c; > 0. 

We then propose and describe the U-CSMA policy as the 
main contribution of this paper. We show that for geometric 
networks, U-CSMA policy overcomes the shortcomings of the 
classical CSMA policy and allows to obtain high throughput or 
utility, arbitrarily close to the optimal, with low packet delay 



that is order-optimal, i.e., stays bounded as the network-size 
increases to infinity. In particular, we analytically show that for 
large networks with torus interference graph and with uniform 
packet arrival rates, the average delay under U-CSMA policy 
is upperbounded as 

Of ' 



independent of the network-size L. 

Using a simulation study, we show that the same general 
delay behaviour also holds for the practical case where 1) the 
arrival rates are determined by a congestion control algorithm 
used on top of the U-CSMA policy to maximize a network- 
wide utility, and 2) the interference graph is geometric (see 
Section |III-A[ ) and constructed in a randomized manner. 

A. Performance of Classical CSMA Policy 

In this section, we provide a motivating example to examine 
the performance of the classical CSMA policy, and explain 
why even for simple topologies, this policy fails to support a 
high throughput with low delay. 

Consider a fixed wireless network with torus interference 
graph, as defined in Section III-C having L links and a 



uniform packet arrival rate A to each link, as defined in d3). 
It is well-known that [38] if all links use the same rate z, 
then the following holds for the achieved uniform throughput 

p,(z,L): 



I h, 



c(L)-v(z > L) = ®(z- 1 ). 



(8) 



This means that to be 0(e) away from the maximum uniform 
throughput /J, m ax(L), an attempt rate z of order - is needed. 

For the above network, two threshold behaviours exist, as 
explained in the following. 

Threshold Behaviour as a Function of Attempt- rate z: It 
is well-known that for a fixed network size L, as the attempt 
rate z increases beyond a threshold, the delay of classi- 
cal CSMA policy on the torus interference graph increases 
substantially. This increase is related to a phase transition 
phenomenon, in terms of the existence of more than one Gibbs 
measures for the infinite torus [39]. 

The currently best explicit characterization of the delay of 
the classical CSMA policy in terms of z shows that the delay 
is (see, e.g., the mixing time analysis in [31]) 

0(z c - L ), (9) 

for some constant c u > 1. While for z < 1, the above 
bound can be moderate for a moderate network size L, for 
z > 1, there will a rapid increase even for moderate values 
of L. Since by ([SI, a large attempt-rate is needed to support 
a high throughput, this explains why the classical CSMA 
policy cannot provide high throughput without incurring a 
large delay. 

We note that by (j8), the classical CSMA policy needs to 
use an attempt rate of order 1/e to support the load p = 1 — e, 
which can be used to write the delay bound in (|9| as 



O 



(10) 



Threshold Behaviour as a Function of Network-size L: 
Depending on the value of a given attempt z, as we increase 
the network size L, the delay of the classical CSMA policy 
shows an undesirable threshold behaviour. 

On one hand, there exists a constant z c .i > such that for 
all attempt-rates z < z c> i, the delay is upperbounded as [42] 



0(log(L)). 



(11) 



This bound states that for low attempt rates resulting in low 
uniform-throughputs, the delay increases only logarithmically 
in the network size L. 

On the other hand, there exists a constant z c ,2 > such for 
any attempt-rate z > z c .2, the delay is lowerbounded as [37] 



ctL/ (log Lf 



(12) 



for some constant q > 0. Hence, for large attempt-rates 
required to support high throughputs, the delay grows expo- 
nentially with the network-size L, which results in a threshold 
behaviour as L increases. It is this exponential increase in 
the delay that prevents the classical CSMA policy to provide 
high throughput with low packet delay as the network-size L 
increases. 

Simulation: To illustrate the threshold behaviours, we have 
simulated a torus of size L 6 {100,400,1600} under the 
classical CSMA policy with uniform attempt rate z. We have 
assumed i.i.d packet arrivals where every unit of time one 
packet arrives for link I, I E C, with probability A, independent 
of any other packet arrival event. For a given network size L, 



to support the uniform arrival rate A (see Section III-C i where 



A = (1 - e)fJ-max{L), 



e > 0, 



(13) 



and consequently a load factor (as defined in |6|) of p = (1 — 
e), we have chosen the attempt rate z such that the resulting 
uniform throughput p(z,L) is given by 



fi(z,L) = p max (L)(l - -) > A. 



(14) 



Fig. |2|a) shows the resulting average queue size per link as 
a function of p in linear scale. This figure clearly illustrates 
the two threshold behaviours. 

First, we see that for a given network-size L, for a small load 
p less than 0.3, the queue-sizes are small. However, as the load 
p increases towards 0.5, which requires a larger attempt-rate z, 
the queue-size increases from only few packets to thousands. 
While the classical CSMA policy is throughput-optimal and 
in principle can support a load p close to 1, we see that in 
practice, it cannot support loads as low as 0.5, i.e., it cannot 
reach the 50% utilization without incurring a large delay. For 
instance, for the 20 x 20 torus, the large delay becomes more 
than lsec for a packet length of 2346 bytes and a channel rate 
of 54Mbs as in 802.11 standards. 

Second, we see that for a given p, the queue-size shows two 
different behaviours. If p < 0.4, the queue-size is small and 
hardly changes with the network size. In contrast, for p > 0.4, 
the queue-size shows a threshold behaviour and drastically and 
exponentially increases with the network size. For instance, at 



p = 0.44, the queue-size almost doubles every time that the 
network size L increases by a factor of 4. 

Intuition: By ([8]), in order to support a high uniform 
throughput, the classical CSMA policy needs to use a large 
attempt rate z. For a large attempt rate z, the network state will 
mainly alternate between two types of transmission patterns 
(valid schedules) where either mostly links in the set of even 
links £( e \ or links in the set of odd links C^°\ are active 
(see Section |HI-C[ ). However, as z and L increase, transitions 
between these two types of patterns occur very infrequently. 
This implies that the classical CSMA policy tends to lock into 
one type of transmission patterns for a very long time before 
it switches to the other type of patterns [39]. 

This locking-in behaviour of the CSMA policy immediately 
implies that while one type of links, e.g., even links, are active 
for a long time, the other type of links, e.g., odd links, cannot 
transmit for a long time. As a result, this locking-in behaviour 
leads to large queue-sizes and hence a large packet delay. 

We next describe U-CSMA policy and provide theoretical 
and simulation results characterizing its performance. 

B. U-CSMA Policy and Its Performance 

The main contribution of this paper is to propose U-CSMA 
policy that overcomes the threshold behaviours faced by the 
classical CSMA policy. As such, U-CSMA enables to obtain 
a high throughput with low delay that is order-optimal. 

U-CSMA Policy: The basic idea behind our proposed U- 
CSMA policy is very simple. U-CSMA policy uses a classical 
CSMA policy z as described in Section III-B However, 
periodically, i.e., at times 

T^iT, ie {0,1, 2,---}, 

U-CSMA policy resets, or unlocks, the transmission pattern of 
the classical CSMA policy by requiring all links to become 
silent, and then immediately restarts the classical CSMA 
protocol to operate as usual. In the rest, we refer to parameter 
T as the unlocking period. We note the in the limit of large 
T, U-CSMA policy reduces to the throughput-optimal classical 
CSMA policy. 

The intuition behind the above unlocking mechanism is to 
prevent the threshold behaviour by preventing the policy from 
locking into a particular transmission pattern for too long. In 
Section [V] we provide one approach to implement U-CSMA 
policy in a fully distributed and asynchronous manner. 

Analytical Results: In order to characterize the performance 
of U-CSMA policy, we first need to know how to choose the 
unlocking period T. While a smaller T helps employ the un- 
locking mechanism more frequently leading to a smaller delay, 
it may also prevent the underlying classical CSMA policy used 
by U-CSMA policy from converging to a maximum schedule 
that is necessary to obtain a high throughput. Hence, as the 
first step, we need to study how fast the classical CSMA policy 
converges to a maximum schedule. 

Our first analytical result (see Proposition [T] in Section VI I 
shows that for the lattice and torus interference graphs with 
uniform attempt rate z, valid schedules under the classical 
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ence graph. ference graphs, as a function of utility ratio p u . 

Fig. 2. Illustration of performance under the classical CSMA policy and U-CSMA policy. 



CSMA policy quickly converge to a maximum schedule at 
a rate that becomes independent of network-size L for large 
networks and attempt-rates. Remarkably, this result shows that 
the distance to the maximum schedules roughly drops as 

1 



Our second analytical result (see Proposition[2]in Sectior VI- 
[B) uses the above convergence result to stabilize networks 
with torus interference topology and uniform packet arrival 
rate A. In particular this result shows that U-CSMA policy 
with unlocking period 



T(e) = 



(15) 



and with large uniform attempt ratd^] z stabilizes the load p = 
(1— e) for large networks with torus interference graph. Hence, 
by the above choice for the unlocking period, U-CSMA policy 
stabilizes queues in the network, all of which have packet 
arrival rate of A = (1 — e)p max (L). 

Further, this result shows that by the above choice for the 
unlocking period T(e), the average queue-size per link and, 
hence, average delay become order-optimal and independent 
of the network size L in the sense that for large L and attempt- 
rate z, they are upperbounded as 



O 



(16) 



Comparing the above delay bound with the ones in (JTOj) 
and ( fT2| i for the classical CSMA policy, we see that U- 
CSMA policy does not suffer from the threshold behaviours. 
Specifically, we see that as a function of 1/e, the queue-size 
under U-CSMA policy increases at most with exponent 3 as 
opposed to the exponent L under classical CSMA policy, as 
suggested by the bound in ( fT0| ). Moreover, U-CSMA policy 
has changed a queue-size that exponentially grows with the 
network size L (see ( fT2| )) to a queue-size that does not depend 
on the network size L. 

4 Large attempt rates can be implemented using Glauber dynamics as in 
[30], [31]. 



Simulation Results: To illustrate the performance of the U- 
CSMA policy and compare it with the analytical results, we 
have simulated a torus of size L € {100,400,1600} under 
the U-CSMA policy. We have assumed i.i.d. packet arrivals 
where every unit of time, one packet arrives to link I, I e C, 
with probability A independent of any other arrival event. We 
have set the uniform attempt rate at z = 50, and for a given 
uniform arrival rate 

A = (1 - e)Hmax{L), 

or load p = 1 — e, we have chosen the unlocking period T as 

1.2 
72"' 



T 



(17) 



Fig. |2|b) shows the resulting average queue-sizes as a 
function of load p. We make the following two observations. 
First, comparing Fig.^b) with Fig. |2ja), we see that while the 
classical CSMA "hits the wall" and its queue-size becomes 
on the order of thousands of packets before reaching a low 
load of p — 0.5, the U-CSMA policy can indeed get much 
closer to the maximum load of 1. In practical terms, for a 
packet length of 2346 bytes and a channel rate of 54Mbs as 
in 802. 1 1 standards, the average packet delay under U-CSMA 
policy becomes 30ms and 90ms for 80% and 85% channel 
utilizations, respectively, while the average delay under the 
classical CSMA becomes more than lsec before even reaching 
the 50% utilization. In addition, replotting the queue-size as 
a function of e = 1 — p in log-log scale (see Fig. Hp, we see 
that the average exponent by which queue-size increases as a 
function of 1/e is 3.02, which closely matches the exponent 
3 as predicted by the analysis in ( [To} . 

Second and as remarkably predicted by the analysis, the 
average queue-size does not change significantly with the 
network size. In fact, for 20 x 20 and 40 x 40 toruses the 
average queue-sizes are hardly distinguishable. This confirms 
that 1) U-CSMA eliminates the threshold behaviours that exist 
for the classical CSMA policy, and 2) the delay under U- 
CSMA is order-optimal in that it stays bounded as the network 
size increases. 

To investigate whether the insight gained through the anal- 
ysis for the torus interference graph carries over to general 



network setups, we have simulated a random geometric in- maintain the length of periods close to T. If t 



(i) 

last 



f (2) 
last 



< T, 



terference graph [2]-[4], see Section III-A in which L G 
{100,400, 1600} links are randomly distributed over a square 
area of 10 x 10, 20 x 20, and 40 x 40, respectively. We 
have chosen the interference range r so that every link on the 
average interferes with six other links. As in [8], [11], [33], 
we have implemented a congestion control algorithm to tune 
the arrival rate to each link so that a network-wide logarithmic 
utility function U ne t is maximized. This algorithm operates on 
top of the U-CSMA policy (see SectionTX|for further details). 

Fig. |2|c) plots the average queue-size as a function of p u 
where 

_ Unet 
Pu — JJ ' 

i.e., p u is the ratio of the achieved network-wide utility to 
the optimum maximal utility U op t- Remarkably, the delay 
behaviour is similar to the one illustrated by Fig. |2|b). 

The main observation here is that the that the insight gained 
through the analysis for the torus interference graph also 
holds for the general case considered here. First, we observe 
that even in random topologies under a congestion control 
algorithm, we can use U-CSMA policy to assign arrival rates 
closely to the optimal without incurring a large delay. For 
instance, for a packet length of 2346 bytes and a channel rate 
of 54Mbs, the delay becomes 40ms to get to 80% of optimality. 
Interestingly, the exponent by which queue-size increases as 
a function of 1/(1 — p u ) approaches 3, the same exponent 
in the delay bound of torus graph in ( fl~6*] > (see Figj7] for the 
corresponding log-log plot). 

Second, we observe that the average queue-size and hence 
the delay slightly change with the network-size. This means 
order-optimality of delay is preserved, and therefore, we can 
use U-CSMA policy jointly with congestion control to assign 
arrival rates close to the optimal with low packet delay in 
arbitrarily large networks. 

In the next section, we provide one example to implement 
U-CSMA policy in a fully distributed and asynchronous man- 



V. Distributed Implementation 

In this section, we provide an algorithm to implement the 
unlocking mechanism of U-CSMA policy, as described in 
Section [IV-B | in a fully distributed and asynchronous manner. 
We assume links can send busy tones to initiate, or relay, the 
unlocking process. Busy tones as opposed to control packets 
are propagated much faster and their detection is easier. 

Fix the unlocking period at T. Let <§; T be the maximum 
delay from the time a link broadcasts one busy tone until all 
its interfering links detect the busy tone. Further, every link 
keeps track of the last time t]^ t and the second last time 

(2) 

Hast tnat e ither initiated or relayed a busy tone. In addition, 
every link maintains a counter that determines the next time 
after irL that it may send a busy tone to locally initiate the 
unlocking process. The value of the counter is reset to T + tf, 
at times wiL}, where is a r.v. Links choose as follows to 



tb is chosen uniformly distributed from [0, 2A], otherwise from 
[— 2A, 0], where A > A;,. A link that joins the network for 
the first time at time t, sets its t^ st = t and its counter value 
to T. 

Every link I implements the following. If at time t, link I 
detects a busy tone or its counter reaches zero, it broadcasts 
one busy tone, to all of its interfering links, only if it has not 
done so in the last 0.5T time-units. This ensures that busy 
tones will not go back to the link where initiated them. After 
broadcasting a busy tone, the link I stops transmitting and can 
start competing for the channel, using the classical CSMA 
policy as usual, only after a time that is uniformly distributed 
in [0, 2A]. This ensures that one link does not always transmit 
first. 

It is clear that for a fixed T and A, as the delay Af, 
approaches zero, i.e., when busy tones propagate very fast, 
the distributed approach converges to the ideal unlocking 
mechanism. For large A;,, however, transmission patterns are 
unlocked locally. Nevertheless, our simulations for both the 
torus and the random geometric networks, as simulated in 



Section IV show that the changes in the queue size and utility 
are less than 1.5% when T = 154, A fc e [0.1, 2], and A = A b , 
all in units of time. Hence, with moderate values of busy tone 
delay, the distributed unlocking mechanism performs close to 
its ideal. 

In the next section, we provide formal statements of the 
analytical results in this paper. 

VI. Performance Analysis 

In this section, we formally state the analytical results de- 
veloped in this paper for lattice and torus interference graphs. 
These results characterize the rate by which the schedule under 
classical CSMA policy converges to maximum schedules, and 
characterize the delay-throughput tradeoff under U-CSMA 
policy. These results use two assumptions that are formally 
stated in Section IVIIII 

Even by making these assumptions, the analysis of CSMA 
convergence is by no means trivial. This analysis requires 
techniques often used to develop mean-field results [40], 
characterizing the properties of ODEs, and also large deviation 



results. Simulation results presented in Section IV and this 
section verify that these assumptions indeed lead to correct 
qualitative results, not only for lattice and torus topologies, but 
also for random geometric networks under congestion control. 

A. Convergence to Maximum Schedules Under Classical 
CSMA Policy 

Our first result characterizes the rate by which the schedule 
under classical CSMA policy converges to maximum sched- 
ules. We consider the lattice or torus interference graph with 
L links, and a classical CSMA policy with uniform attempts 
rate z, as described in Section [HI] 

To state our first result, we use the following notation. Let 
z) be the density, i.e., fraction, of links that are active at 
time t, t > 0. Hence, if N a (t, z) is the total number of links 



that are active at the time t under a classical CSMA policy 
with uniform attempt rate z, then 9i,(t, z) is given by 



;(*>*) 



Ng(t,z) 

L 



We assume that the system is idle at time t = such that 

9 L (0,z) = 0, z>0. 



Let <5k(t, z) be 



5 L (t,z) = 0.5 -9 L (t,z). 



(18) 



Since 0.5 is the fraction of links that can be active under a 
maximum schedule in lattice or torus interference graphs in 
the limit of large L, we see that SL{t,z) can represent the 
distance between the schedule at time t and the limit maximum 
schedules. 

Proposition [T] characterizes how fast the distance 5i,(t, z) 
approaches 0, or in other words, how fast the distance to 
maximum schedules drops to 0, in the limit of large L and z. 

Proposition 1. Suppose the interference graph is given by 
the lattice (or torus) interference graph Gl (or Tl). Under 
Assumptions 7]|2 for Gl (or Tl), there exists a positive 
constant C\, independent of z and L, such that for any r > 0, 
we have that 



lim inf lim inf P 

z— >-oo L— >oo 



sup 

te(o,T 



S L (t,z)- 



Ci 



< 



= 1. 



Proof. Proof is provided in Appendix [TTJ 



□ 



Proposition [T] states that for every finite time-horizon (0, t], 
with probability approaching one as first the network size L 
approaches infinity and then z approaches infinity, the distance 
<5k(i, z) between 6*^(i, z) and the maximum fraction of active 
links 0.5 converges to and drops as O(^) for t £ (0, r], 

The above convergence has two important implications. 
First, under the classical CSMA policy, the distance to maxi- 
mum schedules asymptotically drops as 0(^~), only depend- 
ing on time t. Second, as the O(^) bound does not depend 
on the network-size L or attempt-rate z, the convergence is 
not negatively affected by a large L or large z. This is in a 
stark contrast to the results obtained for for the mixing time of 
CSMA policies, i.e., the rate at which CSMA policies reach 
their steady-state, which increases with attempt-rate z and can 
be exponential in the network size L [37]. 

To illustrate the convergence behaviour, we have simulated a 
nxn lattice, n <E {20, 30, 50, 100}, under the classical CSMA 
policy with z — 100. Fig. [5] shows 9i,(t, z), averaged over 20 
simulation runs, for each lattice. As predicted by Proposition[T| 
convergence behaviour becomes independent of the network 
size for large lattices. In fact, 0z,(t, z) for the largest lattice can 
be very closely fitted by a curve of the form 0. 1(1 + 0. 4i)~ 0,5 , 
which drops to zero as 1/ as stated by the proposition. 
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Fig. 3. Convergence of L (t, z) for lattices of size 20 X 20, 30 X 30, 50 X 50, 
and 100 X 100, from left to right, respectively (z = 100). 



B. Delay-Throughput Trade-off under U-CSMA Policy 

Proposition [T] states that under the classical CSMA policy, 
the distance to maximum schedules converges to zero at a rate 
independent of the network size in the limit of large network 
sizes and attempt rates. Our second result stated in Proposi- 
tion [2] characterizes the delay-throughput trade-off under U- 
CSMA policy for the torus interference graph with uniform 
attempt-rate z (see Section |HI-C| l. Intuitively, Proposition [2] 
states that in large networks, the delay-throughput trade-off 
under U-CSMA policy does not depend on the network size 
L. 

In order to formally state the throughput-delay trade-off for 
any given link in the network, irrespective of its position, we 



consider the torus interference graph 71 (see Section III-Cl 
instead of the lattice interference graph Gl- For the lattice 
interference graph and similar topologies, it is well known that 
due to boundary effects, the throughput achieved by links in 
the network is not uniform over all links in the network when 
a uniform attempt rate z is used [27]. The torus interference 
graph is symmetric with respect to link positions, and as a 
result boundary effects do not exist. While we develop the 
analysis for the torus interference graph, the general insight 
gained through the analysis carries over to more general 
settings, as discussed in Section IV-B 



To state Proposition [2] we introduce several definitions. We 
first note that by Proposition [T] for the torus interference graph 
71 and a given r > 0, we can define a non-negative function 
e p (L, z,t) such that we have 



P 



and 



sup 

*6(0,t] 



S L (t,z) 



< o 



>l-e p {L,z,T), (19) 



lim sup lim sup e p (L, z,t) = 0. 



(20) 



For a given e' > 0, the above limit allows us to define 

z(e',r) and L(z,e',T) such that for z > z(e',r) and L > 
L(z, e', t), we have 

e p (L,z,r) < \d. (21) 
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Fig. 4. log-log plot of average queue-size as a function of distance e to the 
maximum load p = 1, under U-CSMA policy in torus interference graph. 

Furthermore, for a given uniform packet arrival-rate A, < 



A < 0.5, and a given uniform attempt-rate z (see Section III- 
|C) , we define Qi(t, z, A) as the queue size of link Z at time 
i. 

Using the above definitions, Proposition [2] is given as 
follows. 

Proposition 2. Consider the torus interference graph Tl, and 
suppose Assumptions^T^hold for Tl- Let the uniform packet 
arrival rate to each link be A, corresponding to load p{\). Let 
the unlocking period T(X) used by the U-CSMA policy be 



T(A) 



(16d) 2 







where 



e = e(A) = 1 - p(A), 

and C\ is a constant given in Proposition [7] Then, the there 
exists a positive constant C2 such that for z > z(e, T(A)) and 
L > L(z,e,T(X)), the time average of the queue size for any 
link I in Tl satisfies the following under U-CSMA policy with 
the unlocking period T{\): 



lim sup E 

t— >oo 



}l(t,z,X) dt 



C,k^ 



< 



= & 



where k e is defined by the arrival process in Section III-C 



Proof. Proof is provided in Appendix III 



□ 



Proposition [2] states that in order to get e close to the 
maximum load of p = 1, the expected time average of 
any queue-size in the network becomes only 0(fc^/e 3 ), 
independent of network-size L for large L. This is achieved 
by choosing the unlocking period T to be on the order of \. 
By Little's Theorem, we have that the average delay is also 
0(kjL_ /e 3 ). We note that fc c represents the rate by which the 
arrival process converges to it expected value (in the sense of 
Q). Hence, we expect this rate to appear in the average queue 
size and the average delay of any link. In the case where every 
unit of time, packets arrive to each queue according to an i.i.d 
process, we have k e = 1. For such a case, we have that the 
average delay for any given link is 



O 



(22) 
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Quite surprisingly, the above delay-bound and the resulting 
throughput-delay trade-off are valid for arbitrarily large torus 
networks as long as z > z(e,T(X)). Moreover, since C2 
in the proposition is a constant, the delay-bound does not 
depend on the network-size L, and hence, we have an order- 
optimal average delay. This makes the delay-throughput trade- 
off under U-CSMA policy independent of the network-size L 
for large L. As a result, U-CSMA policy, which benefits from 
an unlocking mechanism, can indeed provide high throughput 
with low delay for arbitrarily large torus networks. 

To investigate the accuracy of the delay-bound in ( [22) , we 
have replotted the queue-size as a function of e under the 



simulation setup of Section IV-B The figure shows that the 
queue-size increases with (average) slop 3.02 in log-log scale, 
which, as expected, is close to the exponent 3 given in d22fj 



C. Discussion 

Note that in Proposition [T] and [2] we first let L approach 
infinity and then let z approach infinity. We believe that the 
same result holds if one changes the order of limits. For 



instance, in Section IX we consider a fix network where by 
using different values of the unlocking period, we effectively 
increase the attempt rate. The obtained results, as illustrated 
in Fig. [7] closely match of those if we could change the order 
of limits. We have left a formal proof of this property for the 
future research. 

VII. Dynamics of Schedules Under Classical 
CSMA Policy 

Having provided a formal statement of our main results in 
Section [VI] we now turn our attention to the dynamics of 
schedules under the classical CSMA policy whose characteri- 
zation is the first step for the derivation of Proposition [T] and [2] 

At the heart of the proofs for Proposition [T] and [2] lies the 
analysis of how the density 0L(t,z) of active links evolves 
over time under the classical CSMA policy with uniform 
attempt rate z where all links are idle at time t — 0. To better 
understand the evolution of 0£,(i, z), consider Fig. [5]in which 
9i,(t, z) is plotted for a 100 x 100 lattice. Recalling that each 
unit of time equals to one packet transmission time, we make 
the following observations: 

1) At time t = 0, all links are idle; thus, 9 L (t, z) — 0. 

2) At t = 5, i.e., after five packet-transmission times, since 
the attempt rate z is high, the density z) increases 
quickly to 0.39. 

At time t — 50, the density 6*i(<, z) increases to 0.47. 
At time t = 200, the density f?^(i, z) is slowly reaching 
to the limit of approximately 0.49. 
The evolution of (^(i, z) as explained above is a function 
of the dynamics of CSMA schedules. To see how these 



3) 
4) 



dynamics affect the evolution of z), in Fig. 5(a) Fig. 5(c) 
we have shown three snapshots of the 100 x 100 lattice 

5 As mentioned in Section [iv] we have also observed an exponent close to 
3 when queue-size is plotted against 1 — p u for the case where U-CSMA is 
combined with congestion control in random geometric topologies, implying 
that a variant of Proposition [2] should likely be true for this case. 
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(a) Snapshot at time t = 5, i.e., after five packet-(b) Snapshot at time t = 50, i.e., after fifty packet-(c) Snapshot at time t = 200, i.e. after two hundreds 
transmission times. transmission times. packet-transmission times. 



Fig. 5. Snapshots of the schedules under the classical CSMA policy for the 100 X 100 lattice. A red cross represents an odd active link and a blue circle 
represents an even active link. 



interference graph under the classical CSMA policy at times 
t € {5, 50, 200}. In these snapshots, we have shown only the 
active links. We have shown the even active links by blue 
circles and the odd active links by red crosses (see Section III- 



|C]for the definition of odd and even links). Fig. 5(a) Fig. 5(c) 
illustrate the following typical characteristics of the dynamics 
of schedules under the classical CSMA policy: 

(a) After an initial transient behavior, e.g., at time t = 5 
after five packets transmission times, c/Mifer^] of active 
links have emerged where in each cluster all active links 
are odd and belong to C^°\ or all active links are even 
and belong to C^ e \ 

(b) Shortly after the network starts, e.g., at time t = 5, clus- 
ters of add active links or clusters of even active links 
seem to be uniformly distributed over the interference 
graph. 

(c) Over time, e.g., when we reach time t — 50, clusters 
of odd active links or otherwise clusters of even active 
links start to dominate the schedule. In Fig. 5(b)|5(c) 



clusters of odd active links are dominating the schedule, 
(d) Over time, e.g., at time t = 200, the dominant clusters 
dominate further and grow in size. As a result, the cor- 
responding schedule becomes similar to and approaches 
a maximum schedule in which only one type of links, 
odd or otherwise even, are active. 

Based on the above observations, we can see that Propo- 
sition Q] states how fast the schedule under classical CSMA 
policy consisting of dominant clusters approaches a maximum 
schedule with only odd active links, or only even active 
links, in which half of the links are active. Therefore, to 
prove Proposition [T] as the first step, we need to analyze the 
dynamics of CSMA clusters and their evolution over time. Our 
analysis of the dynamics of CSMA clusters is based on two 
assumptions on the properties of CSMA clusters. We formally 

6 For a formal definition of clusters, their boundaries and areas see Ap- 
pendix [I] 



introduce the two assumptions in the next section. These 
assumptions are made in order to make the analysis tractable 
and we comment on this in more details in Section IVHI-DI 
These assumptions are the only assumptions that we use to 
formally prove Proposition [T] and [2] in Appendix [II] and III 
respectively. 

VIII. Regularity and Randomness Assumptions 

Our analysis leading to Proposition [T] and [2] is based on two 
assumptions on the properties of the clusters that emerge under 
the classical CSMA policy on the lattice or torus interference 
graph with a uniform attempt rate z. The first assumption 
is a regularity assumption which states that the geometry of 
the clusters can not be arbitrary, but satisfy some minimal 
regularity assumption. The second assumptions is a ran- 
domness assumption, which states that while clusters should 
satisfy some minimal regularity assumption, they cannot be 
too regular and need to satisfy some minimal randomness 
assumption. In the following, we first introduce the definitions 
and notations required to formally state the assumptions. 

A. Definitions and Notations 

For the purpose of illustration, we assume that each link I in 
Gl or Tl represented by coordinates (i, j) can be interpreted 
and mapped to the point in M 2 . With such an extension, 
we have mapped the vertex set C of Gl or 71 to a subset of 
points in R 2 . 

We use the following definitions and notations, valid for 
both lattice and torus interference graphs. These definitions 
are formally presented in Appendix [I] For a given cluster C, 
e.g., the cluster of even active links inside the shaded area 
in Fig. [6] we use dC to denote its boundary, and use £(C) to 
refer to the length of the boundary d(C). We also use A(C) to 
refer to the scalar value of the area that cluster C covers. Fig. [6] 
illustrates these definitions. Here, the area that a cluster covers 
and the length of the boundary of a cluster have their usual 
meaning for geometric objects under the euclidean geometry 
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Fig. 6. Demonstration of active links, link coverage area, cluster coverage 
area and its boundary, and bumps in the lattice interference graph G^. Active 
links are colored. 



in M. 2 . We use Cz,(i, z) to denote the set of all clusters that 
exist at time t in the network of size L under the classical 
CSMA policy with attempt rate z, i.e., 

C L (t,z) = {d,l <i< #C L (t,z)}, 

where we define #Ci,{t,z) to be the number of clusters at 
time t. 

At any given time, either the number of even active links is 
the same as the number of odd active links, or the number 
of one type (odd or even) of active links is less than the 
number of the other type (even or odd, respectively) of active 
links. In the first case, we define the non- dominating type to be 
either of the odd or even type. In the second case, we define 
the non-dominant type to be the type of active links whose 
contribution, in terms of the number of active links, to the 
CSMA schedule is less the contribution of the other type. For 
instance, in Fig. |5(c)| even is the non-dominating type. Using 
this definition, we define C^{t, z) to be the set of clusters of 
non-dominating type of active links at time t, and 4t z C i £ d \t 1 z) 
to be the number of such clusters. We define a non-dominating 
cluster to be a cluster that belongs to C^ d \t, z). 

We also define the set Cj™ d \t, z,£) to be the set of all non- 
dominating clusters whose boundary-length is equal to I, i.e., 



C 



(nd). 



L (t, z, I) = {d : d G C ( [ ld) (t, z),t{d) = {,}. 



(23) 



We define #cj" d) (f , z, I) to be the number of non-dominating 
clusters with boundary length t at time t. 

We next state the two assumptions, i.e, the regularity and 
the randomness assumptions. 

B. Regularity Assumption 

The first assumption is a regularity assumption on the 
geometry of a cluster. The intuition behind Assumption [T] is 
that clusters do not prefer a particular direction when they 



grow or shrink, i.e., clusters tend to grow or shrink at similar 
rates in all directions. As a result, it must not be true that 
clusters stay thin so that they have grown only in one direction 
and essentially look like a line. In other words, clusters must 
be fat so that the have grown or shrunk at similar comparable 
rates in all directions. 

Formally, we can define a cluster to be fat if the ratio of its 
area to the square of its boundary length satisfies the following 

A(C) 



e(cy 



(24) 



for some constant c a > 0. For instance, if clusters are 
rectangular and for which, the length of the both sides of the 
rectangles grow by the same factor, as we increase the network 
size L, we then have that the above ratio stays lower-bounded 
at a constant value. On the other hand, if clusters grow in 
only one direction and look like a line, as the network size L 
increases, then the above ratio approaches zero. 

Instead of stating that every cluster is fat so that ( p4| ) 
holds, Assumption[T]states that on the average non-dominating 
clusters are fat so that the ratio 
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is properly lower-bounded, where the ratio is equal to the ratio 
of average area covered by a cluster to the average of the 
square of its boundary length, over all non-dominating clusters 
at time t. Assumption [T] is stated in the following: 

Assumption 1. For the lattice and torus interference graphs, 
there exists a positive constant c a such that clusters under the 
classical CSMA policy with uniform attempt rate z satisfy the 
following: 



lim inf inf 
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The reason why we consider non-dominating clusters is that 
dominating clusters typically are large clusters that often are 
not fat. To see this, consider Fig. |5(b)| in which even (blue) 
clusters are the non-dominating clusters. Intuitively, in this 
figure, even clusters can be considered as "islands" in a "sea" 
of a large odd (red) cluster. While the islands, i.e., even clusters 



in Fig. 5(b) look fat, the "sea", i.e., the large odd cluster in the 
figure essentially contains all cluster boundaries and relative 
to even clusters is thin. 

We also note that in Assumption [T] we first take an inf over 
z and t. This inf makes the ratio lower-bounded independent 
of t and z. We choose to to be bigger than zero since at time 
zero there are no clusters. As far as the analysis is concerned, 
the choice for a positive constant to is arbitrary. Finally, we 
have chosen z > 1 to ensure that clusters exist in the limit of 
large lattices or toruses. 

C. Randomness Assumption 

Assumption [T]requires that clusters (at least non-dominating 
clusters) are regular and fat based on the intuition that clusters 



grow or shrink at a similar rates at all directions. Assumption[2] 
on the other hand requires that clusters are not too regular but 
satisfy some minimal randomness assumption with respect to 
the geometry of their boundaries. To formally state Assump- 
tion [2] we first provide the intuition behind this assumption. 
This intuition does not serve as a proof for the assumption. 

Consider the non-dominating cluster C of even active links 
(inside the shaded area) in Fig. [6] Recall that we use dC to 
denote the boundary of cluster C, as shown by solid lines in 
Fig [6] We see that on the boundary of the cluster, there are a 
number of inactive links. Consider one arbitrary such inactive 
link Iq on the boundary dC, e.g., link Iq in Fig. [6] Consider 
moving along the boundary starting from link l in either of 
the two possible directions, e.g., the one shown in the figure. 
Let U be the ith link visited along the boundary. Suppose after 
visiting n distinct links, we return to link Iq. In such a case, 
we will close the loop and have l n +i = lo- 

We note that moving along the boundary is possible in steps 
of length y/2, as shown in Fig. |6] We define the ith step to 
be the step from link to link l.- L , i > 1, For each step, we 
define its direction to be direction of the movement from link 
link to link ij. For the cluster C, let n r (C) be the number 
steps taken to return to the starting link Iq. Since each step 
is v2 in length, and that the boundary length of cluster C is 
defined as 1(C), we have 



n r (C) 



t(C) 



(26) 



Now, suppose the only information available about a cluster 
C is its boundary length, i.e., we have that 

£(C) = I 

While moving along dC, to close the loop, we inevitably need 
to make direction changes. Formally, we define a direction- 
change event at the ith step to be the event where the direction 
in the ith step is different than the previous (i — l)th step. 
If i = 1, we compare the direction of the first step with 
the direction of the last step returning to Iq. Let pe be the 
probability that there is a direction change clockwise at a 
given step on the boundary. Due to symmetry, pi is also 
the probability of counter-clockwise direction change. Since 
we close the loop, we have that there must be at least four 
direction changes. Hence, for the expected number of direction 
changes, we have 



n r (C)2p e > 4. 
Hence, using ( |26] l, we must have that 



where 



= 2^2. 



(27) 
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(29) 



The probability pg is the probability that there is direction 
change clockwise (or counter-clockwise) at a given step on the 
boundary of a cluster with length £. However, since we close 



the loop, it is clear these direction changes are correlated on 
the boundary of the given cluster. Moreover, these direction 
changes can be also correlated on two different clusters. 
However, as the boundary length I grows, for a fixed i, we 
expect direction changes on the ith step to become independent 
of direction changes at the (i — l)th, (i — 2)th, and the first 
steps. Assumption [2] uses this intuition to assume and state a 
property for the boundary of clusters for all I. 

To state the property, consider moving along the boundary 
dC of the cluster C. While moving along the boundary, we 
may encounter bumps of n steps. Fig. [6] shows bumps of one 
and three steps. Formally, a bump of n steps at a given step 
occurs when starting at the given step, as the first step, 1) there 
is a direction change in the second step, 2) after the direction 
change, there is no direction change in the next n — 1 steps, 
and 3) at the (n + 2)th step there is a direction change so that 
the new direction is opposite to the direction at the first step. 
Hence, when a bump of n steps occurs, for the first time, a 
direction reversal occurs at the (n + 2)th step. 

Now, consider all non-dominating clusters of boundary- 
length I, and suppose an independence assumption holds 
so that a direction change along the boundary of a non- 
dominating cluster with boundary length £ occurs with proba- 
bility pi independently of any other direction on the boundary 
of the same cluster or other clusters. Using this independence 
assumption, the probability of having a bump of n steps at a 
given step is 

Pi- 



Let N^, b \n) be the number of bumps of n steps on the 
boundary of cluster C. Since there are n r (C) steps on the 
boundary, using ( |26] l and ( |28] i, for the expected number of 
bumps of n steps on the boundary of cluster C, we have 



E[7V( fc) (n)] = n r (C) Pl (l - Pl y 
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where 



2c(l- Pl ) 



> 0. 



In particular, for n = 1, we have 



c\ = 2c> 0. 



(30) 
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(32) 



The constant d-y is independent of z, L, or time t. 

As the network size grows, for any fixed £, we expect the 
number of non-dominating clusters with boundary-length £ to 
increase. As such, if a law of large numbers also holds, we 
expect the number of bumps of n steps averaged over all non- 
dominating clusters of boundary-length £ to be at least 
according to ( f30] >. In other words, with probability approaching 
one as L increases, we expect to have that 
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Assumption [2] states that the above inequality holds with 
probability one: 

Assumption 2. For lattice and torus interference graphs, for 
each n > 1, there exists a non-negative constant c n with C\ > 
such that clusters under the classical CSMA policy with 
uniform attempt rate z satisfy the following: 



lim inf inf 
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(34) 



Assumption [2] implies that cluster boundaries exhibit a 
minimal amount of randomness and can not be too regular 
or too smooth. For example, by Assumption |2j it would 
be unlikely to have all non-dominating clusters as perfect 
rectangles, instead the assumption requires that a minimal 
fraction of such clusters of boundary length £ to have bumps 
of for instance length n = 1 along the their boundaries so 
that each such cluster contributes % bumps of unit length 
on the average. As explained earlier, this behaviour would 
be expected if direction changes occurred independently over 
cluster boundaries. This states that the randomness assumption 
can be viewed as a consequence of, and hence weaker than, 
an independence assumption for the direction changes along 
the cluster boundaries. For similar reasons as explained for 
Assumption [T] Assumption [2] takes an inf over z and t. 

D. Discussion 

A few comments on Assumption [T] and [2] are in order. A 
natural question that arises in this context is that whether we 
can prove the conditions of Assumption [T] and [2] for clusters 
under the classical CSMA policy. To answer this question, 
we note that there has been considerable effort in trying to 
derive regularity properties as given in Assumption [T] for 
spatial random processes. But only for processes that are much 
simpler than the CSMA process considered here, such results 
have been obtained. The model for which results have been 
partially obtained, and that is the closest to the random process 
that we consider in this paper, is the Eden model [43] [44]. 

The Eden model is a discrete-time model where initially a 
cluster is given by a single node, and at each time step exactly 
one node on the lattice is added to the cluster boundary, where 
the added node is chosen uniformly and independently (from 
all previous steps) from the set of all nodes that are next to a 
node on the boundary of the current cluster. For this model, it 
has been shown that [43] clusters indeed satisfy the regularity 
condition given by Assumption [T] and we have that 

■A(cy 
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for some constant C. 

The analysis of the Eden model heavily relies on the 
assumptions that 

(a) at each step exactly one node is added to the cluster 
boundary, and 

(b) nodes are added to the boundary of a cluster uniformly 
and independently of each other and all previous steps. 
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These assumptions clearly do not hold for clusters under the 
classical CSMA policy. Lack of these assumptions makes the 
analysis for CSMA clusters difficult. However, it seems true 
that for large clusters, the processes by clusters grow or shrink 
tend become (almost) independent at boundary points that are 
sufficiently far from each other. In such a case, clusters tend 
to grow or shrink in similar comparable rates at all directions, 
which makes it unlikely to have clusters that are too thin. This 
is the intuition behind Assumption [T] 

As mentioned earlier, making Assumption [T] and [2] does not 
make the analysis of the dynamics of CSMA clusters trivial. 
This analysis requires tools and techniques from mean-field 
theory [40], ODE theory [45], [46] , and large deviation results 
[47]. Appendix |U]provides the analysis for the CSMA clusters. 

Simulation results in Section IV and |VI| verify that Assump- 
tion [T] and [2] lead to the correct and precise characterization 
of the classical CSMA behaviour and delay-throughput trade- 
off of U-CSMA policy. As such, we believe that making 
these assumptions is well validated, and the hope is that 
the formulation of these assumptions will also serve as a 
possible starting point for additional studies, possibly allowing 
to further relax or formally prove these assumptions. 

IX. Random Geometric Topologies Under 
Congestion Control Combined with U-CSMA 
Policy 

In this section, we provide the details on how to use U- 
CSMA policy jointly with a congestion control algorithm for 
general topologies. Congestion control is necessary since in 



practice, the capacity region T (see Section IH-B I is often not 
known, and packet arrival rates could initially be outside the 
capacity region T. We provide a detailed look at the simulation 



results provided earlier in Section IV and show that we can 
use U-CSMA policy in arbitrarily large random geometric 
topologies to assign arrival rates close to the optimal utility 
with a low packet delay that exhibits order-optimality in the 
sense that it stays bounded as the network-size increases. 

For simplicity, we assume that links always have data to 
send and consider flows of data instead of discrete-size data 
packets. Using congestion control, we like to ensure that, 1) 
the admitted flows are indeed supportable, and 2) the set of 
admitted flows is chosen in a fully distributed manner such that 
a network-wide utility function is maximized. Suppose Ui(-) 
is the (concave, monotonically increasing, and differentiable) 
utility function for link I as a function of its admitted long- 
term flow rate r\. Suppose the objective is to find \r{\ to have 
the optimal utility U opt [8], [11], [29]: 
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(35) 



Let utility ratio p u be the ratio of the achieved utility U ne t 
to the optimal utility U op t'. 
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Let e u be the distance of p u to the optimal ratio of 1: 

e u = 1 - Pu- 

In order to be 0(y), v > 0, away from U op t, and have e u = 
0(v), one approach [8], [11] is that each link I sets its own 
admitted flow (t) at time i to be 



are max 
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(36) 



where Qz(i) is the queue of admitted flow to link I at time t, 
and £ max is a sufficiently large constant. As for scheduling, 
at any time t, MWM policy can be used that chooses a valid 
schedule I(t) G 1 to solve 



max 



E 



Qi(t)ii(t), 



(37) 



where I(t) = (Ii(t))i e c with = 1 meaning link I is 

active at time f , and (t) — 0, otherwise. The set 1 is the set 
of all valid schedules. 

Since MWM policy is hard to be implemented (see Sec- 
tion [n]), we are interested to use U-CSMA policy for schedul- 
ing. To see how we use U-CSMA policy, first consider a 
classical CSMA policy that sets the attempt-rate of link I E C 
as 



where wi is a weight associated with link I. Let 

W* — max > wiL. 
lex ^— ' 



(38) 



Suppose Proposition [T] extends to random geometric inter- 
ference graphs such that under the classical CSMA policy with 
attempt rates given in ( [38] l and with all links inactive at time 
t = 0, the schedule I(t) used at time t satisfies the following: 



W* 



(39) 



with probability approaching one in the limit of large net- 
works. We note that Proposition [T] can be considered as a 
special case of the above for the torus interference graph 
with wi = 1, I G L. The above extension essentially states 
that we can use CSMA policies to approximate MWM policy 
in random geometric interference graphs, consistent with the 
existence of PTAS for MWM in geometric graphs [2]. 

The above extension motivates us to design U-CSMA policy 
as follows. It resets the scheduling pattern with requiring all 
links to become inactive every T units of time, as described 
in Section |IV-B and sets the attempt rate of any link I at time 
t be zi(t) where 



e kT . 



(40) 



for a fixed integer k. By the above choice, we can ensure that 
the weights do not change substantially^] over an interval of 
length T. At the same time, using the unlocking mechanism 




Fig. 7. log-log plot of average queue-size under a congestion control 
algorithm combined with U-CSMA policy in random geometric interference 
graphs, as a function of the distance e u to the optimal utility ratio p u = 1. 



with period T, we ensure that we never lock into schedules 
for more than T time-units. 

Since we are working with %p instead of Qi(t), we also 
modify the congestion control of ( |36| l so that every link I 
chooses its own admitted flow to be 

Qi(t), 



arg max 

o<a<u< 



2T 



(41) 



Analysis in [8], [11] shows that using the complex MWM 
policy to solve ( [37] > along with the distributed congestion 
control of ( [36} , we will have e u = 0(v) with average packet 
delay of 0(^ _1 ). Using a similar analysis and assuming 
that Proposition [T] can be extended as described earlier, we 
can show that using U-CSMA with attempt rates given in 
( |4"0"| ) and distributed congestion control of ( piT} , we will have 
e u = 0(1/ VT) + 0(f) with average delay of 0{T/v) for 
large random geometric networks. Choosing v = 1/V2T, we 
then have that e u = 0(v) and the average delay as 0(V~ 3 ). 
Hence, to be e u from utility optimality, for large networks, the 
average delay becomes 

CD< <42) 



o 



similar to the delay bound derived from Proposition [2] in ( (22| . 

To investigate the performance of U-CSMA policy, with 
attempt rates given in ( |40] i, used jointly with the congestion 
of ( |4T| , we have conducted simulation for random geometric 
networks of size L G {100, 400, 1600} with interference range 
r such that on the average each link interferes with six other 
links, as described in Section [IV] We have set 

E7i(r) = ]Qg(l + r), I G C. 

We have used different values of unlocking period T, and 
hence, different values of v, in order to obtain different values 
of e u . Note that the exact value of U op t is difficult to compute. 
However, using the fact that log(-) is concave, we can show 
that U op t is upperbounded by Llog(l +Tl), where is the 

7 Attempt rates that are slowly varying functions of links queue-sizes have 
been used in [31], [36] to achieve throughput-optimality for CSMA policies. 



maximum fraction of links that can be activated in a network 
with L links. For the setup considered here, rt approaches | 
for large L. We have used the upperbound for U op t to obtain 
a conservative estimate for e u . 

Fig. [7] replots Fig. |2|c) and shows the average queue-size 
as a function of e u in log-log scale for small e u . We observe 
the following. First, as e u decreases the average queue-size 
increases with a slope close to 3 in log-log scale, as expected 
by ( |42] >, similar to the exponent 3 obtained in Proposition[2]for 
delay-throughput under U-CSMA for torus interference graph. 
This suggests that an extension of Proposition [Tj should likely 
hold. 

Second, we observe that the plots for different network- 
sizes behave similarly. In particular, for large L, i.e., L = 
400 and L = 1600, the average queue-sizes are very close. 
This confirms that U-CSMA exhibits the same delay order- 
optimality and the same desirable delay-throughput behaviour 
observed in the torus interference graph (see Fig. [4]). In 
particular, the simulation results show that we can indeed 
use U-CSMA jointly with congestion control in large random 
geometric networks to operate close to the optimal utility with 
low packet delay. 

X. Conclusion 

In this paper, we have proposed U-CSMA policy as a new 
CSMA policy. In contrast to the scheduling policies in the 
literature, U-CSMA policy not only is simple and distributed 
but also provides high throughput with low delay. Our analysis 
for torus topologies with uniform packet arrivals shows that 
the delay under U-CSMA is order-optimal, and hence, it stays 
bounded as the network size increases. Simulations show that 
the same desirable delay behaviour also holds for the practical 
case where U-CSMA is combined with congestion control 
in large random geometric networks to maximize a network- 
wide utility. Our study in this paper uses a novel approach to 
characterize the performance of random access policies and 
provides a new prospect into the scheduling of large-scale 
multihop wireless networks. 
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Appendix I 

Formal Definition of Clusters, Cluster Area, and 
Length of Cluster Boundary 

In this section, we provide formal definitions of a cluster 
C, its area A(C), and its boundary length 1(C) for the lattice 
interference graph Gl- Corresponding definitions for the torus 
interference graph 71 are defined in a similar manner. 

Given a lattice interference graph Gl, we denote by 8Gl 
the boundary of Gl, i.e. 8Gl is the set of all links for which 
at least one coordinate is equal to or n. 

A. Cluster C 

To define clusters, consider the lattice Gl, and assume that 
links in use the classical CSMA policy with uniform 
attempt rate z for transmission. Recall that in Section III-A a 



link / E C is defined to be active at time t if it is transmitting 
at that time. Then, at any time t, the set of active links can 
be partitioned to a finite number of clusters. We define each 
cluster Cj, 1 < i < i^C^t, z), to have the following properties: 



a) Gj is a subset of active links in Gl- 

b) Connectivity: For any two links I and I', where I ^ V and 
I, I' E Ci, there exists a path of n links {lx, I2, ■ ■ ■ , l n } 
in Ci for some n > 2 where l\ = I and l n = I', such 

that l k = (i k ,jk) E Ci, 1 < k < n, and \i k - i' k+1 \ = 

\h -ifc+il = L 

c) Maximality: Ci is maximal is the sense that no further 
links can be added to Ci without violating one of the 
above properties. 

The above properties define a cluster as a maximal connected 
subset of active links, where each link in the cluster can reach 
any other in the cluster in a sequence of links, or path, in the 
cluster. Considering the mapping from C to M 2 as explained 
in Section VIII-A along the path, the euclidean distance of 
one link to the next is \/2- In Fig. [fi] colored links inside the 
inner polygon represent one cluster. 

By the above definition, each cluster contains only odd 
active links or only even active links. We define an odd (resp. 
even) cluster to be a cluster consisting of odd (resp. even) 
links. 

B. Cluster Area A(C) 

For each cluster C, we use A(C) to denote the scalar value 
of its coverage area Ac- 

To define Ac, first consider links that are inside the lattice, 
i.e., all links / ^ OGl- Any such link has four interfering 
links on the lattice. Considering the mapping from the links 
in C to the points in M. 2 , for any link I inside the lattice, we 
define its coverage-area Ai to be the square-area formed by 
its four closest links. In Fig. [6] we have shown the coverage 
area of one active link. For any link that is not inside the 
lattice, i.e, I E 8Gl, we define Ai to be the intersection of 
the area [0, n] x [0, n] in M 2 and the square that would exist if 
link I were also inside the lattice. For each link I, we define 
Ai to be the scalar value of the its coverage area Ai- 

For each cluster C, we can define its coverage-area as 

A c = U leC Ai, 

i.e., the union of the coverage area of all links that belong to 

C. The coverage-area Ac contains some points where 
links of Gl may be located, and also some points in M 2 where 
links are not located. For instance, the area inside the inner 
polygon in Fig. [6] is the coverage-area of one cluster. 

C. Length of Cluster Boundary £(C) 

For each cluster C, we use £(C) to denote the length of the 
boundary of cluster C, i.e, the length of its boundary dC. 

For each cluster C, we define its boundary dC to be the 
boundary of its coverage-area Ac- Hence, dC is the set of all 
points in M. 2 that any neighbourhood of which contains points 
both in Ac and points not in Ac (see Fig. [6}. 

Appendix II 
Proof of Theorem 1 

In this Appendix, we provide the proof of Theorem 1 for 
the lattice interference graph Gl- The proof for the torus 
interference graph 71 follows similarly. 



Since by ([18), 



S L (t,z) < 0.5 < 1, 



the probability limit in the theorem trivially holds for r < 2 
by choosing C\ > y/2. Therefore, in the rest, we consider 
only the case where 



r>2, 



(43) 



and at the end of the proof, we choose C\ > \/2. 

In the following, we first introduce several definitions that 
will be used throughout the proof. We next in Appendix II 
[A] classify events that occur while the classical CSMA policy 
operates. Using this classification, in Appendix II-B we derive 
a set of stochastic difference equations that characterize how 
the density (fraction) of active links B^t.z) changes over 



time. In Appendix II-C we define the deterministic ODEs 
associated with the obtained stochastic difference equations. 



Finally in Appendix II-D we use several lemmas that are 



provided in Appendix IV to state how the properties of 



the defined ODEs relate to the properties of the obtained 
difference equations, and use these properties to complete the 
proof of the theorem. 

The following are the definitions that will be used through- 
out the analysis. 

Definition 1. The density of events (or links) that satisfy a 
given property is the total number of such events (or links) 
divided by the total number of links L. 

Definition 2. We often define density of links that satisfy a 
given property. To simplify the presentation, we use the defined 
density to also denote the set of links that satisfy the property; 
hence, by writing I £ 6i,(t,z) we mean that link I belongs to 
the set of active links whose density is z). 

Definition 3. For two matrices (or vectors) A and B with the 
same dimensions, we write 

A < B 

if and only if matrix A is component-wise less than or equal 
to matrix B. 

A. Event Classification 

In this section, we classify the events that occur while the 



classical CSMA policy operates, as described in Section III 
[B] This classification provides a basis for the analysis in 
Appendix |II-B| We define four types of events: rare events, 
ordinary events of type-I and type-II, and critical events. We 
start by defining rare events. 

1) Rare Events: We first define the r n -neighbourhood of 
each link I where r n > 10 is a constant. Considering the 
mapping of the vertex set C of Gl to the points in R 2 , we 
define the r n -neighbourhood of any link I to be the set of links 
whose distance from link / is less than r n . 

A rare event occurs at time t if 1) a link / senses the channel 
as idle in the interval [t\,t], t\ < t, and 2) at time t, an active 
link V in the ^-neighbourhood of link I stops transmitting. 



We define links I and V to be involved with the rare event at 
time t. If the links that interfere with link I' find the channel 
as idle when the defined rare event occurs at time t, we also 
define link I and these interfering links to be involved with the 
rare event. Based on these definitions, involved with a rare 
event at time t, there are two or more inactive links within 
r n -neighbourhood of each other that sense the channel as idle 
at time t. We use #L, r (i>z) to denote the density (fraction) 
of links that are inactive at time t, and that in whose rn- 
neighbourhood, there is another inactive link that senses the 
channel as idle at time t such that both inactive links have 
remained idle until time t after a rare event with which both 
links are involved. 

2) Ordinary Events: Depending on the position of links that 
stop transmitting, we consider two types of ordinary events. 
To simplify the presentation, we study these two types through 
the following examples. 

2. a) Casel: Consider the active link a in Fig. [8] that is next 
to a corner on its cluster boundary. We call active links such 
as a that are next to a corner on a cluster boundary as corner 
links. Suppose link a stops transmitting, and suppose no other 
link stops transmitting before link a resumes transmitting 
or link b starts transmitting. In such a case, links a and b 
both sense the channel as idle and compete for transmission 
with attempt-rate z. With rate 2z, one of these links starts 
transmitting before the other. The transmitting link can be 
either a or b with probability 0.5. If link b wins, the corner in 
the boundary moves to north-west, and there is a change in 
the boundary shape. Otherwise, the boundary stays the same. 
A similar discussion holds for links c and d. 

2.b) Case2: Consider the active link e that is next to a corner 
in Fig. [8] We define corners such as the one next to link e in 
Fig. [8] as double corners that are corners on the boundary of a 
cluster extending in two directions at least two steps, each of 
length v2. We define active links such as link e as the double- 
corner links. Suppose link e stops transmitting, and suppose 
no other active link stops transmitting before link e resumes 
transmission or link / starts transmitting. Similar to case 1, 
with rate 2z one of links e or / starts transmitting. If link / 
starts transmitting, which happens with probability 0.5, there 
is a change in the cluster boundary, and two new corners are 
created. 

2.c) Case3: Consider active link g in Fig. [8] Suppose link 



g stops transmitting, and suppose no other link within the 
neighbourhood of link g stops transmitting. In such a case, 
none of interfering links of link g can sense the channel as 
idle. As a result, link g is the only link that senses the channel, 
and with rate z tries to resume its transmission. 

Considering the third case above, we define an ordinary 
event of type-I as the event in which without causing a 
rare event, a non-comer active link, such as link g, stops 
transmitting. Considering the first two cases above, we define 
an ordinary event of type-II to be an event where without 
causing a rare event, a corner or a double-corner link, such as 
link a or e, respectively, stops transmitting. 

We use 9t,\{t,z) to denote the density (fraction) of links 
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Fig. 8. Cases leading to ordinary events. Colored circles represent active 
links. 



that 1) at sometime t-y, t\ < t, have stopped transmitting 
leading to an ordinary event of type-I, and 2) have sensed 
the channel as idle in the closed interval [tx,i\. Similarly, we 
use 0L,2(t, z) to denote the density (fraction) of links that any 
of which, e.g., link a in Fig. [8] 1) at sometime t\, t\<t has 
stopped transmitting leading to an ordinary event of type-II, 
and that 2) the link and the second link associated with the 
ordinary event, e.g., link b in Fig [8] have sensed the channel 
as idle in the closed time interval [ti,t]. 

3) Critical Events: We first define critical events. We then 
explain why these events are critical. 

Consider Fig. [9] in which, we have shown boundaries of two 
neighbouring clusters. Consider the line-segment f — g, from 
link / to link g on the boundary of the right cluster. We see that 
1) this line-segment has unit step length (length of v2), and 2) 
while moving along the boundary, the direction of movement 
before line-segment f—gis opposite to the direction after line- 
segment / — g. Considering this example, we define a critical 
event to be an event in which 1) on the boundary of a cluster 
there is a line-segment of unit step length, and 2) the direction 
before the line-segment is opposite to the direction after the 
line-segment, while the boundary of the cluster is traversed. 
We define an active link such as e as a critical link if around 
which, the boundary satisfies the above two properties. 

To see why critical events are important, consider the critical 
link e in Fig. [9] Consider the following sequence of events, 
which we define as the c- sequence: 

(1) Active link e stops transmitting, 

(2) inactive link / (or g) starts transmitting, and 

(3) inactive link g (or, respectively, /) starts transmitting after 
link / (or, respectively, g) does. 

In this sequence of events, one active link e stops transmitting 
and two new links / and g start transmitting. Hence, the net 
effect is increasing the number of active links by on^] 

Inspired by the c-sequence defined above, we define three 
new densities, all of which as r.v.'s depending on time t. We 
define R,L{t, z) as the density of critical events at time t, i.e., 
the total number of these events divided by L. We use z) 
to also denote the density (fraction) of links that are critical 
at time t, such as e in Fig. [9] 

8 We note that having a reverse sequence of these events, i.e., both links / 
and g stop transmitting and link e start transmitting requires that a rare event 
to occur. 
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Fig. 9. Illustration of the critical event on the boundary of a cluster at 
line-segment f — g. Colored circles represent active links. 



To define 9 L ' (t, z), consider the following. Associated 
with each critical link, e.g., link e in Fig. [9] there are two 
inactive links on the boundary of the cluster to which the 
critical link belongs, e.g., links / and g in Fig. [9] These 
inactive links can sense the channel as idle if the associated 
critical link stops transmitting. We define 9^/^ (t, z) as the 
density of events, i.e., their total number at time t divided 
by L, where both of such inactive links in addition to the 
(previously) critical link sense the channel as idle. Note that 
these events include the event at time t where the following 
has occurred: 1) link e stops transmitting at time t\ < t, 2) 
link / starts transmitting at time t2, t\ < t% < t, and then 3) 
link / stops transmitting at time t so that both links / and g 
sense the channel as idle at time t. We note that with each 
event associated with density 9 L ^ 9 \t, z), three links sense the 
channel as idle. 

Finally, we define 9£ (t,z). As mentioned earlier, associ- 
ated with each critical link, e.g., link e in Fig. [9] there are 
two inactive links on the boundary of the cluster to which the 
critical links belongs, e.g., links / and g in Fig. [9] Consider 
the events where one of these links has become active and the 
other link is inactive. We define 9£ (t, z) as the density, i.e., 
the total number divided by L, of such events. An example of 
such events occurs when the first and second events defined in 
the c-sequence have occurred. In such a case, we have that link 
/ is active and link g is inactive. Since the events associated 
with density 6 L (t,z) involve one active and one inactive link, 
we reuse Ojf (t, z) to denote the density of such active links, 
and also reuse it to denote the density of such inactive links. 

B. Difference Equations for the System Evolution 

In this section, we derive a set of difference equations that 
will be used to characterize the evolution of 9 Lit, z) over time. 
To simplify the presentation, in the rest of the proof, we drop 
dependen cy of the defined link and event densities on z. In 



Appendix II-A we defined densities 0£, r (i), 6i„\(t), 6 l, 2(f), 



R L {t), 9 L m (t), and 9 L f \t). We use these densities to obtain 



the difference equations. 

To obtain the difference equations, we focus on the change 
in the defined densities from time t to time t + e where 

where £ is a constant and < £ < 1. For any fixed z > 0, we 
have that 

lim z 2 e = lim ze — lim zL~^ = 0. (45) 
We obtain the difference equations for 

t>t Q 

where 

to = 1. (46) 

Deriving the difference equations requires us to define r.v.'s 
and obtain their expected values given the history H(t) of 
the defined densities up to and including time t. Formally, 
we define H(t) as the history of all defined densities, i.e., 
L (t), 6 L , r (t), h,x{t), ha{t\ Rdt), 6 ( / g Ht), and 9 L f \t) 
from time zero up to and including time t. 

In the following, we first focus on the change in 0E,(t) from 
time t to time t + e. We then consider the change in other 
densities from time t to time t + e. The analysis in this section 
is based on the following inequality 



o < ( < e L (t) < e u < 0.5, 



(47) 



where 9i and 6 U are positive constants independent of z and 
L. We can assume the above inequality since the theorem is 
stated in the limit of first letting L approaching infinity and 
then letting z approach infinity, and that by Lemma[T]for z > 1 
and all t G [to, r], the above inequality holds with probability 
approaching one as L approaches infinity. 

1) Difference Equation for 9 L(t): In this section, we obtain 
the difference equation for the change in 9 Lit) from time t 
to time t + e given the history H(t). Note that by definition, 
knowing H(t), we know the densities (^(t), 0r, r (t), #L,i(t), 
0~l, 2 (*), R L (t), 9 L fg \t), and 9 L f \t). 

Recall that is the density of active links at time t. 

To obtain the change in 9 Lit), we need to consider all events 
that can affect the number of active links from time t to time 
t + e. The events that contribute to the change in 9 Lit) are 
the following events: 

• A link that is active at time t is not active at time t + e. 

• An inactive link at time t associated with density 9 L ^ (t) 
is active at time t + e. 

• Either of three inactive links at time t associated with 
density 9 L ^ s \t) is active at time t + e. 

• An inactive link at time t associated with density 9L,i{t) 
or 9L,2(t) is active at time t + e. 

• An inactive link I that at time t finds the channel busy is 
active time t+e. This events requires multiple transitions, 
i.e., multiple links changing their states from time t to 
time t + e, in the neighbourhood of the given link I. 
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• An inactive link that at time t is within r n -neighbourhood 
of inactive links associated with density 9L,r(t) is active 
at time t + e. 

In the following, we consider the contribution of each of 
the above events on the change in 9 Lit) from time t to time 
t + e. We start by studying the contribution of active links on 
the change in <?i,(t). 

La) Consider an active link I, and let lg L j(t + e) be the 
indicator function that at time t + e, link I is not transmitting. 
Let Xi(t) denote the rate of change in due to changes 

in states of active links, i.e., let 



Xx(t) 



1 

Le 



le L ,i{t + s). 



(48) 



l<£8 L (t) 



The negative sign is used to indicate that active links that stop 
transmitting decrease the number of active links. We next find 
U[Xi(t)\H(t)]. 

We consider two cases leading to lg L ^(t + e) = 1. In 
the first case, an active link I E 0l(O stops transmitting at 
some time t! , t < t' < t + e, and stays idle by time t + e. 
Recall that packet transmission times and back-off timers are 
governed by memoryless exponential r.v.'s. As a result, since 
links stop transmitting with unit rate, independent of H{t), the 
probability that an active link I G stops transmitting at 

some time t', t < t' < t + e is 



0(e 2 ). 



(49) 



Since the attempt rate is z, once a links stops transmitting at 
time t' , with probability^] 



1 - 0(>£) 



(50) 



it will stay idle until time t + e. Therefore, independent of 
W(t), with probability 



[e - 0(e 2 )] [1 - O(ze)] = e- 0(ze 2 



(51) 



an active link I E stops transmitting at some time t! , 

t < t' < t + e and stays idle until time t + e. 

In the second case, we have lg Lt i(t + e) — 1 as a result 
of link I changing its state more than once. For example, 
this happens when after time t link I stops, starts, and stops 
transmitting all before time t + e and stays idle up to time 
t + e. Independent of H{t), the event that link I changes its 
state more than once, from time t to time t + e, occurs with 
probability not larger than 

[e - 0{e 2 )} [ze - 0({ze) 2 )] = 0(ze 2 ). 

This follows since such an event requires a link to stop and 
start transmitting at least once, each with rate one or z, 
respectively. Hence, the second case occurs with probability 



0(ze 2 ). 



(52) 



9 If a link that has stopped transmitting cannot transmit again, e.g., since a 
neighbour has started transmitting, then the O(ze) term in (50) can be zero. 
The probability 1 — O(ze) accounts for these cases, too. 



(53) 



(54) 



Considering the above two cases, it follows that 

P[le L j{t + e) = 1 \H(t)] =e + e u {t) 

where 

\e u (t)\=0(ze 2 ). 
Having ( |53"| ), we can write 

E[Xt(t) \H(tj\=-j-( Yl ( £ + e i.i(*)) 

iee L (t) 

= -e L (t)(i + ei (t)), 

where 

\e 1 (t)\=0{ze). 
We next find the variance of X\{t) given H(t). Let 
B 1 (t)^X 1 {t)-E[X 1 (t) \U{tj\, 
By ( |53| l and |54|, we have 

P[l 9ll ,(t + e) = 1 \n(t)] < e + 0{ze 2 ). 

Moreover, as shown in Lemma [2j since packet transmission 
times and back-off timers are independent across all links, we 
have that for two different active links I, I' € 0L(t) 



(55) 



(56) 



(57) 



P 



le L ,i(t + e) = 1 |l fli ,i/(* + e) = l,H(t) <e + 0(ze 2 ). 

(58) 



Using the above inequality, |53|, and ( |54] >, we obtain 

^\lg L ^t + e)\e L ,v{t + e)\H(t) 

-E[l 0i ,,(t + e) |W(t)]E[l ei ,j-(i + e) |ft(t)] = 0(ze 3 ). 

(59) 

Using the definitions in (j48j and ( |57| ), we have 



(Lsf 



[Sx(i) 2 |W(t) 

J] E[(l fliii (t + e)-E[l fltll (t + e) |H(i)]) 2 |W(t) 
+ J] E (l flt ,i(* + e) -E[l flj .,,(t + e) 

(l fliiI ,(t + e)-E[l feiI ,(< + e) |W(t)]) |«(t) ^ 

(60) 

Having obtained the above equality, we can use ( |53~| and 
B5J to upper-bound the first and the second summation in the 



above, respectively, leading to 
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E[Bi(t) 2 |ft(t)] 



< 



1 



< 



(Le) 2 ^ 
y ' iee L (t) 

-J_ V 
(Le)* 

1 



(e + 0(ze 2 )) 



0{ze 3 ) 



{LeY 
1 



L^(i)(e + 0(ze 2 )) 
[L<? L (t)] 2 0(z £ 3 ). 



(61) 



Since < 9 Lit) < 1, we can use ( |45[ > in the above to show 
that 

E[Bi(t) 2 |«(*)] < L-t[l + 0{ze)} +0(ze) 

= 7T c + 0(ze). (62) 

Thus far, we have characterized the contribution of active 
links by obtaining ( |55] > and ( |62"1 ). 

7.&) We next consider the contribution of inactive links 
associated with density 9^(t) in the change of 9 Lit) from 
time i to time t + e. Any such inactive link if it is active at 
time t + e, it increases by one the number of active links. 

Recall that 9^\t) represents the density of a subset of 



inactive links (see Appendix II-A.3 1. Consider the inactive link 
I G 9j?\t). For such link I, let ,(t + e) be the indicator 
function that link I is active at time t + e. Using a similar 
analysis used to obtain ( |53"] l and ( |54| ), and recalling that the 
attempt rate is z, we can show that 



P[l gln ,(t + e) = 1 |W(t)] = ze + e 2)J (t), 



where 



|e 2ii (i)| = 0(zV) 



(63) 



(64) 



Define X 2 (i) similar to X\(t) as the rate of change in 9 Lit) 

( f) 

associated with inactive links in 9 1 (t), i.e, let 



Le 



ieei f \t) 



In addition, given let 

B 2 (t)=X 2 (t)--E[X 2 (t) \H(t)]. 
Using ( |63] l and ( |64] i, we obtain 

E[X a (t) \H(t)] = 9 L f \t){z + e 2 {t)), 

where 

|e 2 (t)|=0(z 2 £ ). 



(65) 



(66) 



Using a lemma similar to Lemma [2] and taking similar steps 
leading to d62|, we also obtain 



E 



B 2 it) 2 \U(t) < zL~ c + 0{z 6 e) 



(67) 



l.c) We next consider the contribution of inactive links 
associated with density 9^ it) in the change of 9 Lit) from 



time t to time t + e. By the definition given in Section II-A.3 
9 L ^ (t) represents the density of events where three links are 
inactive and sense the channel as idle at time t. As a result, 
these events each with rate 'dz increase the number of active 
links. This is similar to the previous case where the inactive 
links associated with O/ (t) increase the number of active 
links. Therefore, the contribution of inactive links associated 
with 9 L fg \t) can be characterized in a similar way used to 
characterize the contribution of inactive links associated with 

e { i\t). 

In particular, defining X^(t), similar to X2{t), to be the 
rate of change in 9l^) due to inactive links associated with 
density 9 L fg \t), and letting 

B a (t)=X 3 (t)--E[X 3 (t) \H(t)], 



similar to ( |65| > and ( |66[ ), we obtain that 



E[X 3 (t) \U(t)]=e^(t)(3z + e 3 (t)), 



where 



\e 3 (t)\ = 0(z*e). 
Similar to d67|, we also obtain 



E 



B 3 (t) 2 \H(t) 



< zL 



0(z 3 e). 



(68) 



(69) 



(70) 



l.d) We next consider the contribution of inactive links 
associated with density 9L,i(t) or 0L,2[t) on the change in 
9 Lit) from time t to time t + e. Let 



M = h,xit)+h,2{t)- 



(71) 



Based on the definitions given in Appendix |II- A. 2 the density 



&L,i(t) accounts for ordinary events of type-I leading to only 
one inactive link that tries to access the channel with rate z. In 
contrast, 6^2 (i) accounts for ordinary events of type-II leading 
to two inactive links that with total rate of 2z are trying to 
access the channel. However, in both cases, only one link can 
be made active. As a result, the only difference between the 
contribution of inactive links associated with 0£,,i(t) or 9l,2^) 
is the total rate by which the channel is accessed. In the rest, 
we assume in either case with a lower rate of z the channel is 
accessed. This leads to a lower bound on the fraction of active 
links, which is what is stated in the statement of Theorem 1 . 

Any inactive link associated with 9 Lit) if active at time 
t + e increases the number of active links by one. Define 
X 4 (i) similar to X^it) to be the rate of change in 9 Lit) due 
to inactive links associated with 9l^). In addition, let 

B 4 (i)=X 4 (i)-E[X 4 (f) \H{t)]. 



Taking similar steps used to obtain 
can show that 



where 



E[X 4 (f) \H(t)] =6 L (t){z + tn(t)), 



|e 4 (i)| = 0(z 2 e), 



, (f66l, and d67l), we 



(72) 



(73) 



and 



22 



E 



Bi{t) 2 \U(t) < zL~ c + 0{z i e) 



(74) 



l.e) We next consider the contribution of inactive links 
at time t that sense the channel as busy at time t on the 
change in 9l^), from time t to time t + e. Let 9^ b \t) be the 
fraction of these links at time t. Consider the link I G 9^ it). 
Link I can affect the number of active links only when the 
following occurs. First, all active interfering neighbours of link 
I stop transmitting before time t + e, and second, link I starts 
transmitting before time t + e. For such an event, at least one 



link should stop transmitting, which based on ( 53 1 and ( 54 1 



independent of Wit), occurs with probability at most 

e + 0{ze 2 ) (75) 

Moreover, afterwards, link I should start transmitting, which 
according to (|63"j) and (|64|i, occurs with probability at most 



ze + Oiz 2 e 2 ) 



(76) 



As a result, independent of Hit), the probability that link I £ 
9^ h \t) starts transmitting before or at time t + e is at most 



2_3\ 



ze 2 + 0(V£' 



(77) 



Define X$(t) similar to Xz{Cj to be the rate of change in 
9 l^) due to inactive links that find the channel busy at time 
t. In addition, let 

B B (t)=X B {t)--E[X 6 (t) \H(t)]. 

Treating the density 9i b '(t) in a similar way as we treated 
density 9 L ^\t), and we obtain 

E[X 6 (t) \U(t)] <0<®(t)(ze + e s {t)), (78) 

where 



|e 5 (t)| - Oiz 2 e 2 ). 
Since 9^>\t) < 1, we have 

E[X 5 (i) \H(t)] =0(ze). 
Moreover, we can show that 



E 



B 5 it) 2 \nit) = 0(z 2 e) 



(79) 



(80) 



(81) 



l.f) The final contribution that we consider is due to the 
events that are not previously considered. These events are 
related to rare events, and are the events where an inactive link 
that at time t is within r n -neighbourhood of an inactive link 
associated with density ^L.r(i) is active at time t + e. Define 
X r (t) similar to X\(t) or Xvif) to be the rate of change in 
9 l^) due to these events. In addition, let 

B 6 (t)=X r it)--E[X r (t) \H(t)]. 

Considering the definition for r n -neighbourhood given in 
Appendix II-A.l for any r n , there exists a c„ where c„ 
represents an upperbound on the number of links that are 
within r n -neighbourhood of any given link (. As a result, the 



total number of links that may lead to the final contribution is 
bounded by 



c n L9 L ^ r (t). 



(82) 



Defining 9' L (t) to be the density of links that may lead to 
the final contribution, we have 



9' L M) < c n e Ltr {t). 



(83) 



Treating the density & L (i) the same way we treated density 
9 ( /\t), we obtain 



where 



and 



E[X r (t) \H(t)] <c n e L ^(t){z + e 6 (t)), 



\e 6 (t)\ = 0(z 2 e), 



E 



B 6 (t) 2 \H(t) < zL-t + 0{z 3 e). 



(84) 



(85) 



(86) 



l.g) Having considered the contribution of all events, we 
next derive the difference equation for We first give a 

definition. For any random process h(t) (scalar or vector), we 
define the random process A(h(t)) as 



A(h(t)) 



h(t + e) - hit) 



(87) 



This process measures the rate of change in the random 
process h(t) from time t to time t + e. 

Using the definition of A(-) and the results in |55|l ( 56 1 ([62}, 
(|65j-((6Vj, (^-{TOf, <|72}-(|74]l, (|80)-([8T), (j84}-((86Jl, given 
H(t), we obtain 



A(0lW) 



where 
and 

for which 

E 



- 9 L + 3z9 { L }9) (t) + zO[ n (t) + z9 L (t) 

+ e(t)+B(t) , (88) 

z(t)\ =0(z 2 e) + E[X r (i) \H(tj\, (89) 

6 
i=l 



Bit) 2 \H(t) <0(zL- <: ) + 0(z 3 e). 



(90) 



In the next section, we derive the difference equations for other 
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2) Difference Equations for Other Densities: In the previ- 
ous section, we obtained the difference equation for 6*l(£). We 
can take similar steps to obtain difference equations for other 
densities d^ 9 \t), 9^\t), and 9 L (t). The difference equations 
for these densities are provided in Lemmas [3]|5] 

The difference equations developed in the previous section 
and the ones in Lemmas 3|5 are based on terms such as X r (t), 
Xi f) (t), Xi f9 \t), X r (t), X h (t), and R L (t). In this section, 
we study these terms to obtain a set of difference equations 
that are based on only the densities 9 L {t), 9 {f9 \t), (/) (*)> 
and 9(t). 

Let f2 be the sample probability space and w e O. For a 
given z, we define the event £(L, z, r) as 

£{L, z, r i 





sup 

te[t ,T] 


^L,r{t) < C r Z 2 }n 




{LO 


sup 

te[*o,-r] 


9L,h(i) < c e z~ l }r\ 




{lo 


sup 

te[t ,r] 


9 ( p{t) < cez-'jn 




{lo 


sup 

te[*o,-r] 


9 ( l 9 \t) < C9 z- 1 }n 




{lo 


sup 

te[t ,T] 


9 L (t) < ce« _1 }n 




{lo 


sup 

te[t ,T] 


9 L (t) <0.5-c e , 2 e- 


-}n 


{lo 


inf 

te[t„,T] 


9 L {t) > c 9A e- T }n 




{lo 


ce,i e ~* 


1 < 9 L (t ) < 0.5- 





for some constants c r > 0, eg > 0, c$ t i > 0, and cg.2 > 0, all 
independent of L and z, and for 0L.h(t) defined as the density 
of links that are inactive and sense the channel as idle at time 
t. In the rest, to simplify the presentation, we let 

£{L,z) =£(L,z,t), 

and drop the dependency of £ (L, z, r) on r. Where appropri- 
ate, we do the same for other functions that depend on r. 

Using Lemma [T] Lemma [6] and Lemma [7] we have that for 
any given finite and fixed z > 1 and r > to, 

Urn P\S(L,z)] = 1. (92) 

L— >oo 

Hence, we can define e r (L, z) > such that 

P[£(L,z)]>l-e r (L,z) (93) 
where for any finite and fixed z > 1 and r > to 

lim e r (L,z) = (94) 

In particular, there exists Lq(z) — Lq{z,t) such that for all 

L > L (z) 



e r (L,z) < 0.5. 



(95) 



defined densities 9 ( / g \t), { L f \t), and 9 L (t). 



We next use the definition of £(L,z) to state Ri,{t) as 
a function of 9^(1). Considering ( fl"8"] ) and ( |47] i, we choose 
constants zr and Ljj such that for z > zr and L > Lr if 



£ (L, z) occurs, the conditions of Lemma [I] hold so that w.p.l. 
we have 
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R L (t) > c R 5 L {tf 



(96) 



where c R ; is a positive constant independent of z, L, r, and t. 
Since a larger Ri,(f) indicates more critical events, and hence, 
a larger rate by which 6l (t) increases, in the rest we assume 
w.p.l. 



R L (t) =c R 6 L {tf 



(97) 



Equation ( |97| ) allows us to write RL,{t) as a function of 
(5i(t) and thus as a function of 

We next consider the terms X r (t), xi f9 \t), xi f \t), and 
X r (t) and Xh(t) given in (84) , Lemma[3] Lemma|4] Lem ma [5] 
respectively. We first consider the difference equation in ( |88[ ) 
containing the terms e(t) and B(t), and use the properties 
of event £(L,z) to obtain an upperbound for \e(t)\ given in 
( |89) that contains the expected value of X r (t). To do so, we 
consider the difference equation in ( [88] l conditioned on the 
event £(L,z) as well as history H(t). The equation in ( |88"j ) 
contains some terms such as — and 3z9^ 9 ^ that given 
H(t) are known, and hence, knowing that event £ (L, z) occurs 
does not affect these terms. However, both of e(t) and B{t) 
are affected by knowing that event £(L, z) occurs. 

To upperbound the term |e(t)| given that £(L, z) occurs, we 
first note that by the definition of £(L, z) given in ( |9T) , we 
have 



(98) 



If this inequality was the only information available, in addi- 
tion to H(t), then by ( f84] >, we would have that 

E 



X r (t) |ft(t),gS)J <c n ^ r (t)(z + e6 (t)) = 0(z- 1 ). 

(99) 

However, the event £(L,z) includes other inequalities. To 
address this, we can use the above bound, Lemma [9] and 
assume L > L (z) so that (|95]l holds, to show that 



E 



X r (t) \H(t),£(L,t) 



< 



1 — e r (L, z) 
Thus, given that £(L,z) occurs and using 



^O^ 1 ). (100) 
we have 



\e(t)\ < 0(z 2 e) +E X r (t) \H(t),£(L,z) ^0{z' r ). 

(101) 

For the r.v. B{t) in ( f88] >, using Lemma |9j and assuming 
i > io( z ) so that (|95]l holds, we also have that 



E 



£(t) 2 \H(t),£(L,z) 



E 



< 



-B(i) 2 |W(t) 



1 - e r (L, z) 
< 2[0(zL- c ) + 0(z 3 e)]. (102) 

Hence, 

E\B{t) 2 \H(t),£(L,z)] = 0(zL- c ) + 0(z 3 e). (103) 



We can repeat the same arguments for the difference equa- 



tions for 

<fg)t 



(fa) 



(t), 0£ (t), and 9l(z) containing the terms 



Xy 9 '(t), X}. J >(t), and X r (t) and X h (t), respectively, given 



in Lemmas [3|5] respectively. Doing so and using ( |97| ), and 
applying the definition of £(L,z) to upperbound 0£ (t) in 
3|5 and also to upperbound 0L,h(t) in Lemma [5] 



Lemma 



given 7i(t) and £(L.z), we obtain the following vector 
difference equation for z > z R and L > max(Lg(z), L R ): 



A[9 L (t)] = A9 L (t) + f(0 L {t)) + e(t) + B(t) 
where 



(104) 



e(t) = 





















e(f) 1 


f 


(/«)(*) 




s<»(f) 




e(t) . 



>/(**(*)) 






~c R 8 L {tf 



>B(t) 



flC/)(t) 
. B(t) . 



(105) 



and 



3z 
-3z 
2z 




z 


— z 




z 



— z 



(106) 



Given H(t) and £(L,z), we have that w.p.l 

| e (t)| = 0(z- x ), \eM(t)\ = 0(z- 1 ) 
\eW(t)\=0{z- 1 ), \~e{t)\=0(z- 1 ). 

Moreover, we have 



(107) 



B(t) 2 \%(t),£(L,z) <0{zL-^) + 0{z 3 s), 
BV«\t) 2 \H(t),£(L,z)] <0(zL-<) + 0(z 3 e), 
B^{t) 2 \H(t),£(L,z)] < [0(zL-() + 0(z 3 e), 
B{t) 2 \H{t),£(L,z)] <0{zL- c -) + 0{z 3 e). (108) 



This concludes our analysis to obtain the difference equa- 
tions for the evolution of the defined densities. In the next, we 
define a deterministic ODE as the counterpart of the above 
difference equation. 

C. Deterministic Differential Equations 

In the previous section, we obtained the difference equation 
( |104) for evolution of the vector density 0x(t) from time t 
to time t + e. Using this difference equation, we can obtain a 
deterministic ODE by letting 

||e(t)||=0, ||B(t)||=0, 



and taking the limit of e approach zero while keeping z fixed 



Doing so, and replacing L {t), 9 { / 9 \t), { /\t), and 9 L (t), 



with xi(t), x 2 (t), Xs(t), and x±(t), respectively, and replacing 
5 L (t) with 



6 x (t)= 0.5- an (t), 



we obtain 

d 
dt 

with 



x = Ax+/(x), /(x) 





crSI 




x to = x(t ) = 0l(£q 



(109) 



(HO) 



(HI) 



as the initial condition for x(t) at time io- 

The solution x(i) to the above ODE will be used to 
characterize the evolution of 6l{€) over time. However, we 
need to study the properties of x(t) itself. To do so, we define 



y(t) 



yi(t) 

Va(t) 

V3(t) 
1/4 (*) 



(112) 



to be the solution to the following ODE 

d 



dt 



y = Ay + /(y) 



where 



"0 








" 





-3z 











2z 


—z 





1 








— z 



A = 



with S y (t) defined as 



/(y) = 



l c * S y 


-CRSy 



(113) 



(114) 



S y (t) = 0.5 - Vl (t). 



(115) 



The above ODE has a simpler structure compared to the ODE 
in ( |1 10[ ) since in ( |1 13[ ), given an initial condition, y\{t) can 
be determined independent of y2{t), Uz{t), and 1/4(4). Having 
found the solution for yi(t), one can also find the solutions 
for y 2 (t), y 3 (t), and y 4 (t). 

In order to state how y(t) relates to x(t), we first introduce 
two definitions. We define T>„, as 



= {x:0<Xt< 0.5, i = ,4.}. 



(116) 



Given 2\,, we define a solution x(t) to the ODE in ( |1 10[ ) to 
be well-defined if for an initial condition x fo with x in g P w , 
we also have that 

x(i) G T> w , t G [io,T]. 

In a similar manner, we define y(i) to be well-defined if given 
an initial value for y(t) in T> w , we have that y(i) stays in 



Lemma 10 states how y(t) and x(t) are related, which will 



be used in the final step of the analysis in Appendix II-D 
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D. Final Step 

In this section, we provide the final step of the proof 
for Theorem 1. The goal is to characterize how fast 
diminishes as a function of time t in the time-interval of 
interest (0, t]. Here is the sketch of the final step of the 
proof. We first characterize the initial conditions that will be 
used to define x(i) and y(t) as the solutions to < | 1 1 0| > and 
( |1 1 3| >, respectively. We also define O e =o,B=o(t = ne) as an 
approximation to x(f = ne), where n is a non-negative integer. 
We next show that O^t) is close to y(t) using the following 
inequality 

\\e L {t)-y{t)\\<\\e L {t)-0 L {n t e)\\ 

+ \\0h{n t e) - 0e=o,B=o(n t e)|| 

+ ||0e=o,B=o(rat£) - x(n t e)|| 

+ ||x(n t e) - y(n t e)|| 

+ \\y(n t e)-y(t)\\ (117) 



where for t > 0, 



We properly upperbound each term on the RHS of ( |1 17| i, 
and then use properties of y(t) to complete the proof of the 
theorem. 

To start, we choose L and z such that 



and 



z > max(l, Zr) 



L > max.(L (z),L R ) 



(118) 



(119) 



where Lq(z), Lr, and zr are all defined in the previous 
section. 

Recall that by |92|, for z > 1, the event £(L, z) defined in 
( |9Tj ) occurs with probability approaching one as L approaches 
infinity. The statement in Theorem 1 is in the limit of first 
letting L approach infinity and then letting z approach infinity. 
Hence, in the rest, we assume that event £(L,z) occurs. 
Assuming that event £(L, z) occurs, by ( |9"Tj ), we have that 

< c 6il e- t0 < 9 L (t ) < 0.5 - c e , 2 e- t0 , (120) 

and that for all t € [to, r] 

C£>,ie~ r < 6 L (t) < 0.5 - c e , 2 e~ T , 



and also 



sup max 

teUo.r] 



e[ f9 \t),0^(t),0 L (t) 



< CgZ~ 



(121) 



(122) 



where t = 1 by ( |46] i. 



The inequalities in ( |121[ ) imply that the inequality ( |4"7j ) holds 
which along with ( |1 18[ ) and ( |1 19| l ensures that the difference 
equation of ( |104| > holds for all t E [to, T ]- Moreover, the 
inequalities ( |120| i and ( |122| > at time to and equality ( 111 
provide constraints on the initial condition of the ODE in 
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that define x(t). We will use these constraints to define the 
initial condition for y(t) and to properly upperbound each 
term on the RHS of dTTT*. 



for x(t), by Lemma 10 there exists a t\ independent of z, 
to < t\ < 2 such that 



2h 



Having characterized the constraints on the initial condition Using \\ ( 123| l, ( |126[ ), ( |127| i, (JT28J, and ( |129| l, and 



sup ||x(t) — y(i)|| = 0(z ), as z — » oo, 
tefii.T] 



(123) 



where y(<) is the solution to ( |1 13| > with an initial condition at 
time ti such that 

y(ti)=x(ti). 



By Lemma 10 we also have that regardless of ti, for all 

<G [ti,r] 



y(*) 



and 



0(l)l4xi, as z — > oo, 



(124) 



yi(t)>0.5--^ 

where C y > is a constant independent of ti, r, and z. 
By ( pjg} , and that 



(125) 



It - n t e\ < e, 



we have also have that 



sup || y (rite) - y(i)|| = 0(l)e, as z -> oo. 

te[ti+e,r] 

Define the index n to be such that 

n e {toe~\ toe -1 + 1, io^ 1 + 2, • • • , re" 1 }. 

Define 



(126) 



6> P 



= 0,B 



_o(t = ne) 



to be the deterministic solution to the difference equation of 
( fT04l > with 

||e(t)||=0, ||B(t)||=0. 

Considering the ODE in (JTTOjl, we have that 6> e =o,B=o(^ = 
ne) is a discrete-time approximation of x(t) at times i = ne. 
By Lemma 10 for z > z' Q , where z' Q is a constant, and for t € 
[to, t], we have that x(t) G D w where 2? w is defined in (jTTSJ. 
Therefore, the solution x(t) is bounded over t £ [to,T]. For a 
fixed z > Zq, this boundedness and that /(x) has continuous 
first partial derivatives provide sufficient conditions to have 
(e.g., see Theorem 1.16 in [45]) 

sup ||0 e =o,B=o(n£) — x ( n£ )ll = 0(e), as L — > oo.(127) 

s — — e 

Since we started by assuming that the event £(L, z) occurs, 



by Lemma 11 we also have that for any eg > 



lim inf lim inf 

z— foo L— ¥oo 



p 



sup \\0 L (ne) - O e =o,B=o(ne)\\ < eg 



Finally by Lemma [T2] we have that for z > 1 
lim sup ||0l(*) — ^L(rite)|| = (in prob.) 

L ->°° t£[t ,T] 



= 1. 
(128) 

(129) 



noting that by (|43|> and Lemma 10 



ti < 2 < T, 
we then have that for any r/x > 



lim inf lim inf P 

z— foo L— >oo 



sup ||0 L (*)-y(*)ll <»7i 
te[2,r] 



The above limit implies that 



lim inf lim inf P 

z— »oo L—>oo 



sup |0jr(i) - 2/x C*) | < 7 /i 



1. (130) 



1. (131) 



By ( fT8] l, and recovering the dependencies on z, we have 
<5z(i,z) = S L (t) =0.5-6 L {t). 
By the preceding limit, we then have that 



lim inf lim inf 

z — ^oo L — ^oo 



P 



sup 

*6[2,t] 



6 L (t,z)-[0.5-yi(t)] 



< r/i 



which along with (|125|) implies that 



lim inf lim inf P 

z— foo L— ^oo 



sup 

*e [2,t] 



< r/i 



1, (132) 
1. 

(133) 



Since the choice for r/i > is arbitrary, we can let 

Cy 



which leads to 

lim inf lim inf P 

z— >oc L— yoo 



sup 

te[2,T] 



JL 



8 L (t,z) 



ft 



< 



By defining 



Ci = max(V2,2C v ), 



1. 

(134) 
(135) 



we can extend the above statement for all t G (0, r] since 
always S^(t,z) < 0.5 < 1. This completes the proof of the 
theorem. 

Appendix III 
Proof of Theorem 2 

Here, we provide a drift analysis to prove Theorem 2. By 
the assumption in the theorem, we consider U-CSMA policy 
with the unlocking period T, as described in Section IV-B in 
which transmission patterns are unlocked at times 



Ti=iT, i G {0,1,2,- ••}. 



(136) 



To simplify the presentation, in the rest, we drop the 
dependency of Qi(t,z,L) on z and L. For any link I and 
ti > t\ > 0, its queue size Qi(t) evolves according to the 
following: 

Qi(t 2 ) = Qi(h) - f A(*)lg,(t)>o dt + MhM), (137) 



where 
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AW e {0,1}, 



and Ai(t\,t<z) is the number of packets that arrive to link I 
in the time interval (t\,t%\ as defined in Section III-C| By the 
assumption in the theorem and the definition of e(A) in (IT}, 
the average packet arrival rate to queue Qi(t) is 

A; = A = (1 - e)fi m ax(L) (138) 

where < e(A) < 1. By Q, we also have that w.p.l. 

Mt,t + l) < A 

m ax • 

(139) 

Without loss of generality, for constant A max , we let 
A — 1 

which can be used to show that for t > 0, 

\Ql{t + 1) - Qi(t)\ < m&x{A max , 1) = 1. (140) 

As stated in the theorem, we choose the unlocking period 
T as follows 



T = T(A) = 



(16Ci) ; 

c 2 



(141) 



By ( j!35| >, we have C\ > 1, and hence, for any e, < e < 1, 
we have 



T > 1. 



Let 



(142) 



(143) 



where k t is defined by the convergence property of the arrival 
process for link I in ([3]). Since T > 1 and 

T{i+t)k - Tik = kT > k, 
by the definition of k, we have that 



E 



1 

kf 



-A 



< 



1G 



(144) 



where H s (t) is the system history up to and including time 
Tik- 

We also need to determine z and L. To do so, as stated in 
the theorem statement, we assume that 

z > z(e,T) and L > L(z,e,T), 



where z(e,T) and L(z,e,T) are defined in Section VI-B 
that 



so 



P 



and 



inf (e L (t)-0.5+^=) H 
te[o,T] \ yft 



€ p (L,z,T) < -t. 



> 1 - e p (L,z,T) 



(145) 



(146) 



Having determined the necessary parameters, we next pro- 
ceed with the following drift analysis. Using ( |137[ ) and ( |140) , 
we have that 



Ai(T) = E Qf(T (;l+1) k) - Qi(T ik ) \Qi(T ik ) 
Ai(T lk ,T( l+1 ) k )- 

A(*)lQj(t)>o dt Qi{T lk 



< 2Q z (T lfc )E 



: 



(kT) 2 . 
(147) 



If Qi(T ik ) > kT, by p40] i, it then follows that w.p.l., 

1 Q i (t)>o = 1) * € Pifc,T( i+1 j fe ]. (148) 



For Qi(T lk ) < kT, by ([140}, w.p.l, Qi(T {l+1)k ) < 2kT, and 
hence, 



Ai(T) < 4(fcT) 2 , 
In addition, for Qi(T ik ) < kT, w.p.l., 



i(Ti k ) 



Di(t)l 



Qi(t)>0 



dt < (kT) 2 



(149) 



(150) 



Dividing A^(T) to two expectations, one conditioned on 
Qi(T ik ) < kT and the other on Qi{T ik ) > kT, and using 
the results in ( |147| i-( fT50l l, we can show that 

A<(T) = E[Q?(T (i+1)fc ) - Q?(T lfc ) IQ,^)] 



< 2Q ; (T lfc )E 



z?,(t) |Q,(r< 



5(/cT) 2 . (151) 



We next derive an upperbound for the RHS of (jT3TJ. By 
( |144| >, we have that 

E[MT ik ,T {i+1)k ) |Q,(T ifc )] < fcT[A+^]. (152) 



To obtain a lowerbound for 



E 



T ife 



(153) 



we use the conditions considered for z and L which state that 
([145} and ([146} hold. Since the unlocking mechanism restarts 
the network at time Tj, i > 0, the CSMA transmission events 
after this time are independent of the CSMA transmission 
events before this time. Hence, ( |145[ > and ([146} hold for any 
time interval (7i,Tj + i], i > 0, independent of 7i s (Ti), which 
can be used along with C% > 1 to show that 



E 



1 



> (l-e p (L,z,T)) 



> E 



Ti + 1 



[t)dt H s {Ti) 



0.5 : 



T -1 



T 



2Ci 



> 0.5 - 0.5e p (L, z, T) - [l - e p (i, z)] 



T — 1 

2Ci 



> 0.5 



2Ci 



(154) 



By the choice for T as given in ( |141) , we then have that For t > 0, define 
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E 



(t)dt H s {Ti) 



>0,-|-|-0,-f. 



(155) 

Moreover, since the 2D torus is symmetric with respect to 
link positions, we have that the average transmission, over one 
unlocking period, for any link I is the same as the average 
transmission for any other link I'. This allows us to write 

w e L (t)dt\n s (Ti) 



E 



1 



= E 
_ 1 

~ z 

= E 



lec 
Tj Tl 



L 



^Di(t)dt H s (Ti) 



T 



EC 
T i+1 



Di(t)dt HsiT,) 



D l {t)dt % s {Ti) 



(156) 



for any I € C. By the above equality and ( |155| >, we obtain 



E 



T 



D t (t)dt \H.(Ti) 



Ti 



> 0.5 



(157) 



Using ([138]), < [I5T) , ( fT52| ), ([157), and that n max {L) < 0.5, 
we have that 



\(T) < 2Q{T lk ) \kT{\ + ie) - fcT(0.5 - |e) 
L 16 8 

= -e'kTQ(T lk ) + 5(kT) 2 



where 



1 



-6 < 1. 



5(fcT) 2 
(158) 

(159) 



Summing over i = {0, • • • , J — 1}, and then taking the 
expected value, we obtain 



E 



i=0 



= E\Qf(T Jk )-Qi(T = 0) 
j -i 



< -e' 



'E[^ferQ z (T ife )] + J5(fcT) 2 , (160) 

i=0 

Rearranging terms in ( |160| ), dividing by JkTe', letting J 
approach infinity, and assuming queues are initially bounded, 
we obtain 

,7-1 



1 / \ 5A-T 

< I™ sup (E [Q 2 (0)] + J5(fcT) 2 ) = — . 

(161) 

From (1140b and that fc > 1 and T > 1, we obtain 



r ( i + i) fc 



)l{T ik ) + {\t\-T ik ) 



dt 



< kTQ l {T lk )+2{kTf 



(162) 



it 



r * 



By the non-negativity of queue sizes and that t < Ti tk , we 
have that 

1 /"* ukT 1 [ Ti * k 

- J o Qt(t) dt< * t . tkT j o Qi(t) dt, (163) 

which along with (|162[) leads to 



1 



Qlit) dt < 



i t kT 1 
t itkT 



u-l 



kTQi(T mk )j +i t 2{kT) 2 
(164) 

Using ( | 1 6 lb , taking the expected value of both sides of ( |164[ ) 
and then the limsup as t approaches infinity, and using the 
limit 



lim 

t->oo t 



lim 

t— ¥00 



t 



we obtain 



lim sup E 

t— >oo 



Qi(t) dt 



< 5- 



kT 



Using ( |141| i and ( |159| l, we then have that 

56fcT 



lim sup E 



Qi{t) dt 



< 



1. 



2kT. 



C 2 k 

e 3 



(165) 



(166) 



where the constant C*2 is given by 

C 2 = 56(16Ci) 2 . 



Recalling that fc is defined in ( |143| l, we have completed the 
proof of Theorem [2] 

Appendix IV 
Lemmas 

Lemma 1. Consider r > to and suppose z > 1. There exist 
positive constants cg i and eg 2 , independent of t, z, and L, 
such that 



lim P 

L— too 



< inf 6 L (t) 
te[t ,T] 



n( sup 6 L (t) < 0.5 - ce i2 e T 



Moreover, for t > 



lim P 

L— f oo 



ce,ie 
Proof. Recall that 



to 



<0 L (t o ) < 0.5 - c e ,2e- t0 



t a = 1. 



1. (167) 



= 1. (168) 



Consider link I, which has a maximum of four interfering 
neighbours in the set Mi. Recall that at time t — 0, all links 
are idle. Let time wi be theirs? time after time t = that link 
I or one of its interfering links in Mi starts transmitting. 

We claim that at any time t < wi, at least one of the links 
in the set {1} UMi senses the channel as idle, independent of 
states of other links in Gl from time to time r. To reach 
a contradiction, suppose t < Wi, and it is true that all links 



in {1} U Mi sense the channel as busy at time i. Since the 
first transmission occurs at time wi, this means that link Z 
is not transmitting at time t and finds one of its neighbours 
transmitting at that time. A neighbour transmitting at time f 
means that by definition wi < t, contradicting the assumption 
that t < wi. 

Define the set as the set of links in the (n+l)x(n+l) 
lattice Gl with coordinates of the form 



(4i,4j), < 4i < n, < 4j < n. 



The set is a subset of even links C^ e \ By the definition link I is not transmitting, as defined in Appendix II-B.l and 
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for some constant eg 2 > independent of r, z, and L. 
Combining the last two limits, we obtain the first limit in the 
lemma, as required. The second limit also follows by replacing 
t with to in the above discussion, completing the proof. □ 

Lemma 2. Consider two active links I and I' at time t, i.e., 
suppose 1,1' £ 0l{£). We have 

le L ,l(t + e) = l \l eLtV (t + e) = l,H(t) <e + 0(ze 2 ) 

where lg L j(t + e) is the indicator function that at time t + s, 



of £.i e \ we have that for any two links 1,1' € Cf , the sets 
{1} UMi and {/'} [J Mv are disjoint. 

By definition, at time toj, for the first time link I or one of 
its interfering links starts transmitting. Let t; the amount of 
time that such a first transmission lasts. Since before time wi, 
at least one link in the set {1} U Mi is trying to access the 
channel with rate z, and since packet transmission times are 
independent of each other with unit rate, we have that 



P 



< %ti>r \{w v ,t v ,l' ±1,1' e£W} 



U'l 



>(1 



-0.5z 



)e-\ 



(169) 

Let I T j be the indicator function that the for the first time 
after time zero link Z or one of its interfering links starts 
transmitting at time Wi < 0.5to = 0.5, and that such a 
transmission continues until and including time r. For z > 1, 
by (|169|) we have that 



P 



/r,i = i \{U,i',l' ±1,1' eC^} 



> Pr 



(170) 



where 



p T = (1 - e-°- 5 )e- T 

Therefore, we have a sequence {I T .i, I € ci^} of indicator 
functions that independently of each other will be one with 
probability at least p T . Similar to the W.L.L.N, it then follows 
that 



> 0. 



lim P 

L—too 



E 



led 



Ir.i 



> 0.5p T 



1. 



(171) 



By definition, 6 Lit) is the fraction of all active links at time 
t, and for t 6 [ioi T ]> it includes all links I 6 df* 1 for which 
I T j = 1. Hence, for t € [to >■?"]> w.p.l 



;(*)> 



Letting 



C0,i =0.25(1 -e" - 5 ), 
and using the above and (|171[), we obtain 



lim P inf e L (t) > c e ie~ T = 1. 

L->oo LiG[t Q ,r] 



Using similar arguments, we also have that 



lim P 

L— > oo 



sup ^i(t) < 0.5 — ce, 2 e T 
*e[t ,T] 



= 1, 



(172) 



(173) 



"H(t) w f/?e history of densities as defined in Appendix 
Proof. By ( |53"j ) and ((54]), we have 

P[l fl£ ,i(t + e) = 1 |W(*)] - e| = 0(ze 2 ). (174) 
We also have that 



P 



l fllll (t + e) = 1 and U L ,i'(t + e) = 1 \U{t) 



< P link Z and l' stop transmitting before t + e \ H(t) 
= [e - 0(e)][e - 0(e)] <e 2 . (175) 
The first inequality follows from the fact that having 

±e h ,i{t + e) = 1, 

by definition, requires link Z to stop transmitting before time 
t + e. The equality follows form ( |4"9| ) and the fact that packet 
transmission times are i.i.d over time and over links. 

Using the definition of conditional probabilities, ( |174[ ) and 
( |175| l, we obtain 



P 



< 



,/(t + e) = 1 \lg L .y{t + e) = 1 \H(t) 
e 2 

n , 2 , =£ + 0(ze 2 ), 
- U{ze z ) 



as required. 



□ 



Lemma 3. For the density 9^ 9 ^ (t), we have that 
A[9[ f9 \t)} = \Rl (t) - 3z^ /9) (t) + 9 {} \t) 



(176) 



where 



\e if9 \t)\ = 0(z 2 e) + I][X^ 9 \t) \H(t)], 
E[X^)(t) |W(t)] <c„0 iir (t)[z + O(z 2 e)]. 



(177) 



anc/ c n > M a constant independent of L and z as defined 
in Appendix \II-B.l\ f. 

Proof. The events that can contribute to the change in O^ 9 ^ (t) 
from time t to time t + e are the following events: 

• A critical link at time t associated with density i?i(t) is 
inactive at time t + e. 



• An active link at time t associated with density 6 L ^\t) 
is inactive at time t + e. 

• Either of three links at time t associated with density 
9/ (t) is active at time t. 

• Multiple transitions (state-changes) from time t to time 
t + e done by one link or a link and the links in its r n - 
neighbourhood that can increase or decrease the number 
of events associated with 9 L ^ 9 \t + e). 

• An inactive link that at time t is within r n -neighbourhood 
of another inactive link associated with density 0L,r(t) is 
active at time t + e. 

Considering the above events and taking similar steps leading 
to |[88|l-(|90), we obtain the statement of the lemma, as required. 

□ 

Lemma 4. For the density 9 L (t), we have that 
&[9[ f \t)} = \2z9[ f9 \t) - (1 + z)0W(t) 



+ e^\t) + B^\t) 



(178) 



where 



\e^(t)\ = 0(z 2 e)+nx ( r f \t) \H{t)\, 
EpT r W(t) \U{t)] < c n 9 L At)[z + 0(z 2 e)] 
E[B^(t) 2 \H(t)] <0(zL-() + 0(z 3 e). 



(179) 



and c n > is a constant independent of L and z as defined 
in Appendix \II-B.l\ f. 

Proof. The events that can contribute to the change in 9^ (t) 
from time t to time t + e are the following events: 

• Either of the two of the three inactive links at time t 
associated with the events defined by 9 L ^ 9 \t), e.g., link 
/ or g in Fig. [9] is active at time t + e. 

• An active link at time t associated with density 6 L (t) 
is inactive at time t + e. 

• An inactive link at time t associated with density 0£ (t) 
is active at time t + e. 

• Multiple transitions (state-changes) from time t to time 
t + e done by one link or a link and the links in its r n - 
neighbourhood that can increase or decrease the number 
of events associated with 6 L *\t + e). 

• An inactive link that at time t is within r„ -neighbourhood 
of another inactive link associated with density ^i, r (t) is 
active at time t + e. 

Considering the above events and taking similar steps lead- 
ing to (|88|l-(|9"0"l), we obtain the statement of the lemma, as 
required. □ 

Lemma 5. For the density 9l (t), we have that 

A[9 L (t)} = [(9 L (t) - R L (t) - 9^{t)) - z~9 L {t) 

+ e(t)+B(t)], (180) 



where 

\e{t)\ < 0{z 2 e) + nx r (t) \H(t)]+E[X h (t) \H{t)], 
E[X r (t) \H(t)} < c n 9 Wr (t)[z + 0(z 2 e)], 
nx h (t) \U(t)]<c n 9 L . h {t)[l 
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0{ze)}, 



E 



B{t) 2 \H(t) <0{zL-t) + 0{z 6 e), 



(181) 



and 9L,h(f) is defined in Appendix II-B.2 and c n is a constant 



independent of L and z as defined in Appendix II-B.l f. 



Proof. Recall that Otif) (see Appendix 



pendix II-B.l 



II-A.2 



and Ap- 



represents the density of links that are inactive 
at time t as a result of an ordinary event of type-I and type-II. 
Considering this, the following are the events that contribute 
to the change in O^t) from time t to time t + e: 

• An active link I at time t is inactive at time t + e where 
link I at time t 1) is not associated with densities or 
ffv'(t), and 2) at time t is not within ^-neighbourhood 
of any inactive link that senses the channel as idle at time 
t. 

• An inactive link at time t associated with density #z,(i) 
is active at time t + e. 

• Multiple transitions (state-changes) from time t to time 
t + e done by one link or a link and the links in its r n - 
neighbourhood that can increase or decrease the number 
of events associated with 9 Lit + e). 

• An inactive link that at time t is within r n -neighbourhood 
of another inactive link associated with density 9L :T (t) is 
active at time t + e. 

^L,h{t) is the fraction of 



As defined in Appendix 



II-B.2 



links that are inactive and sense the channel as idle at time t. 
Considering this definition, we note the following. In the first 
event considered above, we included active links not associated 
with densities Rhif) and 9^\t) since these active links are 
already accounted for in the difference equations for 9^ 9 '(t) 
and 9^\t), respectively. We further included those of such 
active links not within ^--neighbourhood of inactive links 
that sense the channel as idle to avoid rare event^] By the 
definition of c„, the density of links within ^--neighbourhood 
of inactive links that sense the channel as idle is not larger 
than 

c n 6 L ,h{t)- 

However, if an active link / is within ^-neighbourhood of an 
inactive link /' that senses the channel as idle at time t, and 
1) link /' does sense the channel as busy in the time interval 

[t' ,t + e], for some t', 

t <t' <t + e, 
and 2) link I stops transmitting at some time t" , 

t' <t" <t + e, 

10 We need to avoid rare events since the definitions of densities 0L,l(t) 
and (>L 2(t) requires active links stop transmitting without leading to rare 
events. 



and stays idle by time t + e, then we may have an ordinary 
event. These sequence of events are accounted for by multiple 
transitions event considered above. 

Considering the above observations and taking similar steps 
leading to (f88]>-([90|), we obtain the statement in the lemma, as 
required. □ 

Lemma 6. Consider any finite titne-interval [ii,^]. where 
< <i < ti < oo. Then, for z > 0, 



lim P 

L— >oo 



sup 

te[*i,*2 



,r{t) < 



1 



where c r > is a constant independent of z and L. 



Proof. Recall that in Section II-A.l we defined #£ r (£) to 



denote the density (fraction) of links that are inactive at time t, 
and that in whose r n -neighbourhood, there is another inactive 
link that senses the channels as idle at time t such that both 
inactive links have remained idle until time t after a rare event 
with which both links are involved. 

In addition, recall that $l,/i(£) is the density of links that 
are inactive and sense the channel as idle at time t. By ( |194[ ) 
in the proof of Lemma [7] we have that 

A{9 L<h {t)) <[5-{z + 5)9 L>h (t)]e + e h (t) + B h (t) (182) 



where e/j(£) and Bh(t) are properly bounded in ( |192| i and 
( |193[ ), respectively. 

To obtain a difference inequality for 0L iT -(£), we study the 
change in #L )r (£) from time t to time t + e. Consider a link 
I G 0L,h(t)- If an active link /' within 2r„-neighbourhood 
of link I stops transmitting at some time t' > t, then at 
most link I, link I', and interfering links of link I', which 
may also find the channel as idle, may be counted as the 
links belonging to 0L,r(t')- We chose the 2r„-neighbourhood 
instead of r n -neighbourhood towards obtaining an upperbound 
for the rate by which $L,r(£) increases. Therefore, since links 
stops transmitting with unit rate, and that there are at most c' n 
links in the 2r„-neighbourhood of any link, where c' n > is 
a constant, the term 

Qc' n 9 L , h (t) 

plus an error term gives the maximum rate by which 0L,r{t) 
increases due to active links at time t that are within 2r„- 
neighbourhood of inactive links in 0L,h(t) an d that are idle 
by time t + e. 

In addition, since each inactive link / in 0L,r(t) with rate z 
tries to become active, the term 

z6 L , r (t) 

plus an error term gives the minimum rate by which 9^ r (t) 
decreases due to inactive links in 0L,r(t) start transmitting. 
Finally, we must consider the events in which multiple tran- 
sitions occur from time t to time t + e by one link or a link 
and the links within its r n -neighbourhood. These events can 
contribute to 9z r (t + e). 
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Considering the maximum and minimum rates obtained 
above and the events with multiple transitions, we can use 
a similar approach as taken in Appendix |II-B| to show that 

A[9 L , r (t)} < Sc'JlA*) - z0 L At) + er(i) + B r {t). (183) 
where 



\e r (t)\=0(z 2 e), 



and 



^[B r {t) 2 \H{t)\ = 0(zL~<) + 0{z A e) 



(184) 



(185) 



For a fixed z, we can obtain the deterministic differential 
inequality counterparts of ( |182| > and ( |183| l by letting 

e n (t)=B n (t) = e r (t) = B r (t) = 0, 



letting 



0. 



and replacing 0L,h(t) and 0L. r {t) with x^it) and x r (t), 
respectively: 



dt 
d_ 

dt 



Xh < 5 - (z + 5)xh, 



(186) 



By the proof of Lemma |7J we can choose t! x independent 
of z and L with < t\ < t% such that 



sup x h (t) < c' g h z \ 

*6[*i,*2] 



(187) 



where c' g h > is a constant independent of z and L. Using 
this upperbound for Xh(t) and the differential inequality for 
x r (t) in ( |186| l, we can show that 



sup x r (t) < c' r z 2 
te[*i,t 2 ] 



(188) 



for some constant c' r > independent of L and z. 

Therefore, the solution to the deterministic differential in- 
equality counterparts of ( |183| l is less than c' r z~ 2 over the entire 
interval [£i , £2]- Using this result, the bounds given in ( |184[ ) 
and (|185[), which both approach zero as L approaches infinity, 



and the methods of Appendix II-C and II-D we obtain the 
statement of the lemma for 



1. 



as required. 



□ 



Lemma 7. Consider any finite time-interval [£i,£a], where 
< t\ < £2 < 00. Then, for z > 0, 



lim P 



sup 

te[ti,t 2 



max (0 Ljh {t), 9 { L f9) (t), { L f) (t), 9 L (tj) < cgz^ 1 
where eg > is a constant independent of z and L. 



Proof. We first provide the proof sketch for the statement of 
lemma only containing 6L.h(t), i-e., we show that 



lim P 



sup 0L.h(t) < ce,h,z 
te[*i,t 2 ] 



= 1 



where eg /, > is a constant independent of z and L. 
Taking similar steps, we obtain the same limit probabilities 



for 6j[ (t), 6V'(t), and 0E,(t) with corresponding constants 



3 (/) 



> 0, 



> 0, and eg > 0. Since all of these limit 



probabilities state that events occur with probability one in 
the limit, the intersection of these events also occurs with 
probability one in the same limit as L approaches infinity. 
Hence, using these limit probabilities, we can obtain the limit 
probability in the lemma with 

c g = max (ce,h, 4 , c g f l ,Cg). 



We now proceed to provide the proof-sketch for ( |189) , By 
definition, since 0L,h(t) is the fraction of links that are inactive 
and sense the channel as idle, and 0l (t) is the fraction of links 
that are active, both at time t, we have 



(190) 



Consider the change in 0L,h(t) from time t to time t + e. 
Since the attempt rate is z, using the same arguments provided 



in Appendix II-B we can show that the term 



z9 L ,h(t)e 

plus an error term gives the fraction of links that belong to 
0L,h(t) and that are active at time t + e. Hence, z9l h,(i)e plus 
the error term serves as a lowerbouncp] for the rate by which 
&L,h(t) decreases from time t to time t + e due to inactive 
links in &L,h(t) that start transmitting. 
Similarly, the term 

59 L (t)e 

plus an error term gives the maximum increase in QL,h{t) due 
to active links in Oi,{t) that stop transmitting. This follows 
since 1) each active link stops transmitting with unit rate, and 
2) at most five links may sense the channel as idle when an 
active link stops transmitting, i.e., the link itself and its four 
interfering links. Finally, we must consider the events in which 
multiple transitions occur from time t to time t + e by one 
link or a link and the links within its r„-neighbourhood. These 
events can contribute to LJl (t + e). 

Considering the obtained rates and the events with multiple 
transitions that contribute to the change in O^j^t) from time 
t to time t + e, we can use a similar approach as taken in 
Appendix |II-B| to show that 

A(0£,j,(i)) < -z0 L ,h(t)e + 59 L {t)e + e h (t) + B h (t) (191) 

"We have a lowerbound since an inactive link in 0L.h{t) mat ' s active at 
time t + e may make other (interfering) inactive links in 9^ ^ (t) sense the 
channel as busy at time t + e. In such a case, by defintion, QL,h{t + e) w iU 
not account for these inactive links. 



where 
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\e h {t)\=0{z 2 e) 1 



(192) 



(189) and 



nB h {t) 2 \H h {t)} = 0(zL-<) + 0(z 3 e) (193) 

where Hh(t) is the history of 0L,h(t) from time zero up to 
and including time t. Using ( |190| > and ( |191| l, we have that 

A(0z, h (i)) < [5 - (z + 5)9 LJl (t)]e + e h (t) + B h {t) (194) 

For a fixed z, we can obtain the deterministic differential 
inequality counterpart of the above stochastic difference in- 
equality by letting 

e h (t)=B h (t) = 0, 

taking the limit of e approaching zero, and replacing 0£ /,(£) 
with Xh(t): 

d 

—x h < 5 - (z + 5)x h . 

Any solution to the above inequality is bounded by a function 
of the form 

-^-+C h e^\ 
z + 5 

for some constant Ch > that depends only on the initial 
condition on ^(O). Since as the initial condition 



we have 



< x h (0) = 6 L<h (0) < 1, 



< C h < 1. 



Therefore, considering the assumption in the lemma that < 
ti < ti < oo, we can find a constant c' g h > independent of 
L and z such that 

5 



sup Xh(t) < sup 
te[ti,t 2 ] te[ti,t a ] 



-(z+5)t 



<c 



8,h* 



(195) 



Therefore, the solution to the deterministic differential in- 
equality counterpart of ( |194[ ) is less than c' g h z~ Y over the 
entire interval [ti,^]- Using this result, the bounds given in 
([192} and ( |193| l, which both approach zero as L approaches 
infinity, and the methods of Appendix |II-C| and |ITD[ we then 
can show that for 

c e,h = c' gh + 1 

( |189| > holds, as required. □ 

Lemma 8. Consider the lattice interference graph Gl and a 
time t, to < t < oo. Suppose we have that 

< 5i < S L {t) <S 2 < 0.5, (196) 
for some constants Si and 62- Moreover, suppose 

rj(L)+59 LJl (t) <0.5S L {t) (197) 
where rj(L) is defined in Lemma [77] with the property that 
lim rj(L) = 0, 



and 0L,h(t) is the density of links that are inactive and sense 
the channel as idle at time t, as defined in Appendix \II-B.2\ 
Then, there exist constants z' R and L' R such that for z > z' R 
and L > L' R , w.p.l, we have 

R L (t) > c R 5 L {t)\ (198) 

for some constant cr > 0, independent of t, L, z, and r. 

Proof. By the assumption in the lemma, we have that 

r){L) + S6 L , h {t) <0.5(S L (t). (199) 

Using this inequality, Lemma [To] and the definitions given in 
Section [Villi we have w.p.l 

^^cec\ 



> 1. 



(200) 



L6 L (t) 

Let £(t) be the average boundary-length of non-dominating 
clusters at time t, i.e., let 



£{t) 



E 



cec 



(201) 



Using this definition and ( |200| i, we have that 

#C L nd \t,z)£(t)> -j=L8 L (t). (202) 

Since the total number of links in Gl is L, and the coverage 
area Ai of each link is at most two (see Appendix |TJ|, we find 
2L as an upperbound for the total area covered by the links 
in the clusters. Hence, we have that 



E 

cec l L nd \t,z) 



A(C) < 2L. 



(203) 



By Assumption [TJ for any z > 1 and t > to, there exists Lai 
such that if L > Lai, w.p.l, we have that 



]T A(C) > 0.5c o 



cec£" d) (M) 



cec 



E 



£{Cf 



(204) 



(M) 



In the rest, we assume that 

L > Lai, 

and 

z > z' R = 1. 
By ( p03| ) and ( pCM) , for L > L A i> we obtain 



E 



4L 



cec^ a '(t,z) 

By Jenson's inequality and ( |201| >, we also have that 



E 



cec 



...-) (t „,/(C)' 



#4 nrf) (*^) 



> 



E 



cec 



which along with ( |205| ) leads to 

#C { [ ld \t,z)£(tf< 



AL 



(205) 

£(i) 2 , 
(206) 

(207) 



Using \202) and \2Q1) , we have that 

4^2 1 
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£{t)< 



Using the above inequality and (|202|i, we have that 



(208) 



(209) 



We use the above bounds to find a lowerbound for i?£,(t). 
Recall that R,L(t) is the density of critical events (see Ap- 
pendix II-A.3[ ). Since each critical event maps to a bump of 
one step (see Appendix VIII-Ci, is also the density of 

bumps of one step. By Assumption [2] w.p.l, for a given £, for 
non-dominating clusters of length I, we have that the lim inf 
of average number of bumps of one step per non-dominating 
cluster is lower-bounded by ^h. Therefore, for L > Li, where 
Lg is a constant dependent on £, we have that 



E cecr ) (t , z , ^ h) (- = i) 



> 



0.5ci 



(210) 



>(t,z,£) * 

where c\ is a positive constant independent of z, t, r, or L. 

Consider all non-dominating clusters of boundary length 
k£ (f) or less where 

k = 2. 



By \\96) and ( |208| l, £{t) is finite and bounded. Hence, there 
are only a finite number of values for the boundary length 
of the considered clusters. This implies that we can find a 
constant Lai such that if L > La2, then ( |210| i holds for all 
£ < k£(t). In the rest, we assume that 

L > L' R = max(L j4 i, L A 2) 



so that ( p08| ), ( p09l >, and <j2T0j> hold for I < k£{t). 

By ( |208| l and Markov inequality, for the total number of 
non-dominating clusters with boundary length k£(t) or less, 
we have 

]T #c£*>{t,z,t) >(l~l) *C ( [ ld \t,z). (211) 

l<kl(t) 

Using ( |208> , p09l , ( |2T0] >, and (|2TTj, for the total number of 
bumps of one step on the boundary of non-dominating clusters 
with boundary-length of k£(t) or less, we have 

E E ^ h) (D> E #4" d) (^,^ 

^ e #cr(M^) 

v ' t<kl{t) 

>^ r (l-\)L5 L {tf. 
~ kUV2\ kJ v ' 



(212) 



Thus, for 2 > z^, and L > L' R , we have 



R L [t)>\ E ^'W^AW 3 (213) 



t<kt(t)cec£ d) (t, z ,i) 



where 



is a constant independent of L, z, r, and t. This completes 
the proof of the lemma. □ 

Lemma 9. Consider a non-negative r.v. x, and an event A 
such that P(A) > 1 — ea, where < < 1. We have 



E[x\A] < 



l-e A 

Proof. Define A c to be the complement of the event A. For 
the expected value of r.v. x, we can write 

E[x\ = E[x\A]P(A) + F,[x\A c }[l - P(A)}. 

Using the non-negativity of r.v. x and the assumption for the 
event A, we have 



E[x\A] 



Efxl 



E[x\A c }[l - P(A)} < E[x] 



P(A) 



P(A) 



< 



completing the proof. 



□ 



Lemma 10. Suppose the following holds for the initial con- 
dition x to = x(fg) of the ODE of (JTTOjl. First, suppose 



< c Xl<1 < xi(t ) < c Xu 2 < 0.5, 



(214) 



where c Xlt \ and c Xlt 2 are constants independent of z. Second, 
suppose for 2 < i < A, 



x l {t ) = O{z~ 1 ). 



(215) 



Then, there exists a constant z' such that the following holds 
for z > z' . First, we have that the ODE of < | 1 1 0| > has a well- 
defined unique solution x(i), i.e., we have that x(i) exists, is 
unique, and x(i) G T> w for all t G [ioi T ]- Second, there exists 
t\, to < t\ < 2, independent of z, such that 



sup \\x(t) — y(t)\\ = 0(z 1 ), as z 
te[t!,T] 



(216) 



where y(t) is the solution to ( |1 1 3| > with the initial condition 

y(<i)=x(i 1 ) 

with the property that 
d 



sup 
and for t>t\ 



df 



y{t) = 0(l)l 4x i, as z 



yi(t)>0.5-^L 



oo, 



(217) 



(218) 



where C y > is a constant independent of t\ and z. 

Proof. By Lemma [13] there exists a constant zq such that for 
z > zo, there exists a ii independent of z with 

h < h < 2 
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such that under the assumptions in the lemma, x(i) is unique 
and well-defined over the interval [to-^iL i- e -> 

x(t) €V W , te [t ,h}. (219) 

Lemma [13] also states that for z > zo, we have 

0<4 ljl <x 1 (t 1 )<4 1>2 <0.5, (220) 

where d x 1 and c' 2 are constants independent of z. By 
Lemma [T3| we also have that 



Xi{h) = 0(z- 1 ), 2<i<4, 



and 



3z 
2 

1 



c R 5Uh 



0(z- 2 ) 
0(z~ 2 ) 



x A {h) ~ - \xtih) - cuS^h) 6 = 0(z- z ) 



(221) 

(222) 
(223) 
(224) 



Having that for z > z$, x(i) uniquely exists over the interval 
[tojti], it remains to prove that x(t) as the solution to 



—x = Ax + /(x), x(fi) = x tl 



(225) 



is also unique and well-defined over the interval [ti , r] given 
an initial value at time t\, and that ( |216[ )- pT8| hold. We note 
that general theorems on existence and uniqueness of solutions 
for ODEs require Lipschitz continuous functions on the RHS 
of ODEs. However, the RHS of \225\ is only locally Lipschtiz 
continuous. One approach to address the lack of Lipschitz 
continuity is through adding a term to the RHS of ( |225| l to 
make sure that the RHS is Lipschtiz continuous. To do so, we 
define the following ODE with an initial condition at time t\, 



where t\ is specified in Lemma 13 

d 



dt x a = Ax 3 + /( x a) + S(x g ), 
x s (*i) = x(£i) 



(226) 



where x s (£) is given by 



Xg,\{t) 
Xg,2{t) 

x g,3(t) 

Xg,i{t) 



and we choose g(x) to be a function that satisfies the follow- 
ing: 

. For ||x|| > 2, 

||/(x)+. 9 (x)|| = 0. (227) 

• Defining 

T> = {x;xi e [-0.6,0.6],l<i<4}, (228) 
we have that 

|| S (x)||=0, xeX> (229) 



The function (/ + <?)(x) is continuously differentiable 
(on M 4 ) and Lipschitz continuous with Lipschitz constant 
Lf +g that is independent of z. 

The Lipschitz continuous function (/+<?) (x) is such that 
for any xi and X2 

+ (/ + .g)(x 2 ) = L /+S (xi,x 2 )(xi -x 2 ) 

where each element [Lf +s (xx,X2)]jj of matrix 
Ly +s (x 1 ,X2) is such that 



|[L /+9 (xi,x 2 )]y| < L 



f+g- 



(230) 



The third condition in the above ensures that for a given 
initial condition at time t\, the solution x g (t) to the ODE 
in ( |226[ ) exists for all t > to and that it is unique (e.g. see 
[46], Chapter 3, Theorem 3.). We also note that a function 
g(x) satisfying the fourth condition exists since /(x) is a 
(polynomial) function of only x\. 

Consider the ODE for y(t) as given in ( |1 13| > with the initial 
condition at time t\ such that 



y(tx) = ^(ti) = x(i x ). 



(231) 



Considering the constraints on the initial condition y(£i) as 
imposed by the above equality and ( |220| i-( |224] i, it follows from 
Lemma Q3] that 



sup 



dt 



y(t) 



o(i)u 



(232) 



By the ODE for y(t), we also have that yi(t) is given by 
the following ODE 



c R 6*{t) = -c R (0.5- yi ) 3 



dt a 3 
with the initial condition that 

Vi(ti) = xx(ti). 

Using this ODE for yi(t), we obtain for t>t\ 

y 1 {t)=0.5-a(l + p(t-t 1 ))-i, 

where 

a = 0.5 - yi (ti), 

and 

P 

By ( |220T > and ( p34) , we have that 

< a% < a < a.2 < 0.5, 
< ft < /3 < ft < Or, 



4c fl 2 



(233) 



(234) 



(235) 



(236) 



(237) 



(238) 



where ai, a 2 , ft, and ft are constants independent of z. Using 
these constants, we can find a constant C y > independent 
of ti and z such that 



yi (t) > 0.5 



(239) 



We can also rewrite the ODE for y\(t) given in ( |233| l as 
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2/1 = -yi + 3zy 2 + zy 3 + zy 4 

-3zy 2 + c R 5l(t) 
2 

- ^2/4 +2/1 - c R S^(t). 



(240) 
(241) 
(242) 



By Lemma 15 we have that the sup of expressions in ( |240| i- 
( |242| > over [ti,r] all are 0(z _1 ). Hence, for t € [*i,t], we 
have 



dt 



Hi 



-yi + zy 3 + zy 4 + e y 



where 



sup \e yi (t)\=0(z- 1 ). 
te[*i,T] 



(243) 



(244) 



Since by Lemma 15 for t € [ii,T] we have y(t) £ T> w , 
and O w C T>, by the property of g(-) given in ( |229| l, we have 
that 

llff(y(*))ll=0, *G[*i,r]. 



Using this equality, the ODE in ( |243| l, and the differential 
equations for y2(t), yz{t), and y4(t), as given in ( |1 13[ )-( fTT4"| >, 
we obtain that for t G r] 



y = Ay + /(y)+ S (y) + e y , 



(245) 



where e y (t) accounts for the error term e Vl {t), and by ( |244| i, 
we have 



sup ||e y (t)|| =0(z- 1 ). 

te[ii,r] 



(246) 



Therefore, y(t) is the solution to a perturbed version of the 
ODE in ( |226| l. We next use this result to complete the proof 
of the lemma. 
Define e g (t) as 



e g (t)= Xg (t)-y(t). 



(247) 



Considering the properties given for the chosen function g(-) 
in the beginning of the proof, by Lemma [14] we have that 



sup ||e 9 (t)|| ^Oiz- 1 ). 
te[t!,T] 



(248) 



Since by Lemma 15 y(t) £ T> w where V w is strictly inside 
T), by ( |248| l, we have that for z > z\, where Z\ > z and is a 
sufficiently large constant, 



Xg(t)€V, te [h,T] 



Hence, for z > z x , by ( g29J ) and p49| >, 

g( Xg {t)) = 0, t G [tx,r\, 
which along with (|226]) for z > Z\ leads to 

J^g = AXg + f(Xg), t G , t] . 



(249) 



(250) 



(251) 



In addition, since yi(t) is an increasing function of t and 
less than 0.5 as given by ( |235| l, by ( |248) >, we can choose a 
sufficiently large constant z<z, z-i > z\, such that for z > z-i 
and any t E [ti,r], 



-(0.5 -a 2 ) < x g ,i(t) < 0.5 



(252) 



where «2 is a constant defined in ( |238| l. Using the above 
bounds for x g ,i(t) and that x ff (i) satisfies the ODE in ( |251[ ) 
with the constraints at time t\ given in ( |220| i-( |22"Tj ), we have 
that 



< x g ,i{t), 2 < i < 4, te [ti,r] 



and 



sup a; fl ,i(f) = 0(z x ) 
te[*i,T] 



2 < i < 4. 



(253) 



(254) 



Therefore, by the above bound and ( |252| i, for z > z 3 where 
Z3 > Z2 is a sufficiently large constant, we have 



x fl (t) eD w t e [fi,r] 



(255) 



Since x g (i) is uniquely defined over the interval [ti,r], 
by ((226), j25l\ , and ( |255| >, we have that for z > z 3 , x 9 (t) 
serves as a well-defined solution for the original ODE in ( |225| >. 
Moreover, since the RHS of the ODE in ( |225| l is continuously 
differentiable with respect to Xj's, 1 < i < 4, we are ensured 
if a solution to ( |225| l exists, it has to be unique (e.g., see 
Theorem 1.13 in [45]). Therefore, x g (i) serves as a unique 
solution to ( |225| l for the interval [ti,r]. By ( |248| l and ( |255] >, 
it then follows that for 

z > z' Q = z 3 , 

the ODE in ( |225| l has a well-defined unique solution over the 
interval [ti,r], and that 



sup ||x(i)-y(<)|| =0(z~ 1 ), 



(256) 



which along with the results in ( |232| i and ( |239| l completes the 
proof of the lemma. 

□ 

Lemma 11. Given that event £{L, z) as defined in ( |91[ ) occurs, 
for any eg > 0, we have that 



liminf liminf 

z— voo L— >oo 



P 



sup \\0 L {ne) 



0e=o,B=o(ne)|| < eg 



Proof. Consider the system at discrete times in the set 

{ne, n > toe^ 1 }. 

For these times, by definition, ^(rte) provides a solution to 
the difference equation in ( |104| >, and e =o,B=o(«£) is the 
solution to the same difference equation with 

||e(t = ne)\\ = 0, ||B(t = ne)|| = 0, 



and the initial condition that 

^e=0,B=o(io) = ^i(io)- 



Mi 



(257) 



Throughout this proof, by the assumption in the lemma, we 
assume that event £(L, z) occurs. As a result, all expectations 
and probabilities are stated conditioned on £ (L, z). In the rest, 
to simplify the notation, where appropriate we write n to 
indicate the time instant t = ne; hence, we write 6l(ti) to 
mean 9i(t = ne). 

We first provide the proof-sketch. Define the sequence 
{T n ,n > < £ -1 } as 

n—toe -1 

T n — £ I |e(n — m)\ + |B(n — m) 

where e(n) — e(t = ne) and B(n) = B(t = ne) are defined 
in ( fTUB) . Define 



=0,B 



--o(n). 



(258) 



e e (n) = L (n) - 
By ((257), we have 

e g (n = toe^ 1 ) = 0. 
In the first part of the proof, we show that if 

sup T„ < ei 

where t\ > can be chosen arbitrarily small, and if 

||0 e =o,B=o(i = ne) - x(< = ne)\\ < eg, — < n < - (261) 

e e 



(259) 



(260) 



where x(i) is the solution to ( | 1 10[ >-( |TTT| > and eg > can be 
chosen arbitrarily small, then for z > z' Q , where z' Q is defined 
in Lemma fTOl we have 



sup ||e e (n)|| < eg. 



(262) 



In the second part of the proof, to complete the proof, we 
show that ( |260| l holds with probability approaching one as we 
first let L approach infinity and then let z appraoch infinity. 
We also show that \26\\ holds for sufficiently large L. 

First part: Suppose ( |260| > and \26\\ hold for arbitrarily 
small e\ > and eg > 0. We choose eg arbitrarily small 
such that < eg < 1, and let 

1 



(\ = £0 



4ce ic 



(263) 



where c is a constant independent of z that will be determined 
later (in the following inductive proof). We next show that 

to 



IMn)|| <4ce L ^T„, 
which along with ( |260[ ) and ( |263| > leads to 



< n < -. 

e e 



sup ||e e (n)|| < eg. 



(264) 



(265) 



What remains to complete the first part is to show that 
(|264)> indeed holds. We first note the following. Subtracting 



the difference equation associated with <?(ri) e =o.B=o from that 
of 01, (n) given in ( |104| i and using ( |258| l, we obtain 

A[e e (n)} = Ae e (n) + f(8(n)) - /(0(rc) e=o ,B=o) 
+ e(n) +B(n), 

If || eg (n) || < 1 and ||0i(n)|| < 2, we can write 

A[e e (n)] = Ae g (n) + L(n)e e (n) + e(n) + B(n), (266) 

where L(n) is a 4 x 4 matrix such that 

|L y (n)|<i C) 1 < i,j < 4, (267) 

for some constant L c > 0. This follows from the definition of 
/(•) as a polynomial function, as defined in ( |105[ ). Therefore, 
as long as 0(n) and eg(n) are properly bounded, the evolution 
of eg(n) can be locally linearized using the matrix L(ri) as in 
fl266l >. 

We next prove inductively that for all 

n e {toe^ 1 < n < re" 1 }, 



1) inequality ( j264| > holds, 2) equality ( |266| ) holds, and 3) 
0(n) e V where V is defined in ( |228] >. 

For k = trjE -1 , by ( |259| l, we have the degenerate case of 

e e (fc) = 

for which ( |264[ ) holds. Moreover, since the event £(L,z) 
occurs, we have that 9^(1 = to) € T>, or in other words, 
9h{k = io£ ) € X>. This implies that 

||0 L (A : = t o£ - 1 )|| <2. 

Since e e (fc = t £ _1 ) = and ||0l(A: = io^ 1 )!! < 2, we 
also have that equality ( |266| l holds. Hence, for fc = toe _1 , the 
claims of the induction hold. 
Next, suppose for some fc where 



^o^ 1 < k < re" 



1. 



the claims of the induction hold for all 

n € {toe^ 1 ^qe^ 1 + 1, • • • , fc}. 

We show that the claims hold for time n = k + 1. By the 
induction assumption, ( |266| l holds for all 

n G {i e _1 ,ioe _1 + 1 ■ ■ • 

For such an n, using the definition of A(-) given in ( (87] i, we 
can rewrite ( |266| > as 

e e (n + 1) = A(n)e e (n) + e[e(n) + B(n)] (268) 

where 



Using the above and ( |259| l, we have thaj*^] 
e e {n + l) 

n—to£~ 1 m— 1 

- E [II M 
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m=0 



J=0 



Using d267b, we then have that 



n— toe 

|e fl (n+l)|<e £ A' r ' 

m— 



e(rt — m) + B(n — to) 

(270) 

|e(n — m)| + |B(n — to)| 

(271) 



where 



A = I + e[A + £ c l 4 x4]- 



(272) 



All eigenvalues of A are distinct and real, and for small 
e, they are close to one. Moreover, we have that the largest 
eigenvalue of A is less than 1 4- 4eL c . Using this and the 
eigenvalue decomposition for A to write it as A = UDU -1 , 
where D is a diagonal matrix containing the eigenvalues of 
A and U is the matrix with its columns as the eigenvectors 
of A, we can show that each element A™* of the matrix A™ 
satisfies the following 



| A™ | <c(l + AEL c ) m <ce 4L ^ £m , 



(273) 



for some constant c > 0. This constant was used earlier to 
determine e± in ( |263| l. 

Using the bound in ( |273| l, we have that 



|e(n — to) | + |B(n — to)| 
<4ce L " Em |e(n-TO)| + |B(n-TO) 
Hence, by ( |271| i and ( |274| l, we have that 

||e fl (n+l)|| < 4ce i ^ n+1 )T„. 
Since, by the induction assumption, 

n e {t s~ 1 ,t e~ 1 + !,•■• ,k}, 



(274) 



(275) 



we can use choose n — k in ( |275[ ). Therefore, we obtain that 
p64) holds for n = k + 1. 

Moreover, by ( |260| ), ( |263| >, and the inequality ( |275| ) for n = 
k, we have that 



\e e (k + 1)11 < e e < 1. 



(276) 



In addition, by ( |261| l and the error bound in ( |276| l along with 
definition of e$(n) in ( |258| >, we have that 

\\0 L (t = (k+ l)e) - x(i - (fc + l)e)|| < 2e e . 

Under the assumption in the lemma that £(L,z) occurs, by 



Lemma 10 for z > z' , we have 

x(t = (n+ l)e) e A, 



A(n) = I + e[A + L(n)]. 



(269) 



l2 For m = 0, we define the product term in (270} to be one. 



where T> w is strictly inside T>. This and the above inequality 
imply that for small eg, we have 

9 L (k + l) G V. 

This inequality implies that ||0(A; + 1)|| < 2, which along with 
( |276[ ) confirms that ( |266| l holds for n = k + 1. Therefore, we 
have that for n = fc + 1 all claims of the induction hold, 
completing the inductive proof. 

Second part: By the non-negativity of the norm function, 
the sequence {T„} is a submartingle. By Doob's inequality, 
for any €\ > 0, we then have that 



P 



sup T n >e 1 \£(L,z) < ej^EfTj \£ (L, z)] . (277) 



By ( |107| i, we have that 

||e(n)|| =0(^ 1 )- (278) 

In addition, by ( |108[ ), the following inequality for two given 
random vectors a and b: 

E[||a+b||] <E[||a||] +E[||b||], 

and that for a given r.v. x, E[x] < yTEja^, we have that 



E 



||B(n)|| \£(L,z) < O(z - 5 J L-°- 5<: )+O(z 1 ' 5 e - 5 ) 



(279) 



From ( j278| ) and ( |279| >, we have that 

r - 1 



E 



Thus, 



\£{L,z) 



1 



O(z- 1 )+O(z - 5 L- a5C ) 



0(z 



1.5 £ 0. 5) 



(280) 



limsup limsup E[Ti |£(Z, jj)] = 0. 



L— >-oo 



Taking the limit of both sides of ( |277| >, we obtain that for any 

ei > 



lim sup lim sup P sup T„ > t\ \£(L,z) 



= 0, (281) 



which means that 
lim inf lim inf P 

z— yoo L— too 



sup T„ < ei \£(L, z) 



= 1. 



(282) 



The above limit states that ( |260| l holds with probability ap- 
proaching one in the limit as first L approaches infinity and 
then z approaches infinity, as required. 

We next show that ( |261| l holds for sufficiently large L, By 
the discussion in Appendix |II-D| leading to ( |127) , having that 
the event £(L, z) occurs, we can show that 

sup ||0 e =o,B=o(ne) — x ( n£ )ll = 0(e), as L — > oo. 

(283) 

Therefore , for any eg > 0, we can choose L sufficiently large 
so that ( |261| l holds, completing the second step and the proof 
of the lemma. 
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Lemma 12. Suppose z > 1. For t > 0, define 

t 



n t 



We have that 



lim sup \\e L (t)-e L (n t e)\\ =0 (inprob.) (284) 

L ^°°te[to,r] 

Proof. Consider the evolution of the system in the time 
interval [ne, (n + where n G {0, 1, 2, • • • }. Suppose for 
a given n, there exists some t G [ne, (n + l)e) such that 



-6> L (ne)|| > 7 /2 , 



(285) 



where 772 > 0. 

Since is a four dimensional vector, the above inequal- 

ity implies that within the time interval [ne, t], at least one of 
the densities in the vector O^t) must change by an amount of 
Therefore, a fraction C772, where c > is some constant, 
of links must change their states at least once, from active 
to inactive or vice versa, from time ne to time t. Define 
Xi(n) £ {0,1} as follows. Let xi(n) = 1 if link I changes 
its state at least once from time ne to time (n + l)e, and we 
let Xi(n) = 0, otherwise. 

Based on the above discussion, (|285[> implies that 



tec 



xi(n) > cLn 2 



(286) 



We note that a^'s may be correlated random variables. To 
address this difficulty, we focus on the times epochs after time 
ne that one link changes it state. Since packet transmission 
times and back-off periods are continuous r.v.'s, w.p. 1 ., no two 
links can change their state at the same time. This implies that 
the amount of time from time ne until the first link changes 
its state is well-defined, which we denote it by si. Hence, 
the first state change occurs at time ne + s%. Similarly, let Sj, 
i > 2, denote the amount time from the (t — l)th state change 
to the 7th state change. Hence, the 7th state change after time 
ne occurs at time 



ne - 



3=1 



Without loss of generality in the rest we assume that cLr\i 
is an integer. Assuming so, the condition in ( |286| l implies that 
the cLrj2th state change after time ne must occur up to time 
(n + l)e, and therefore 



cLrj2 



Si < e. 



(287) 



Consider the time ne + 2Jj=cT s j- At this time, there are 
a number of links that sense the channel as idle and try to 
access the channel with rate z, and there are a number of 
links that are busy with packet transmission and with unit 
rate stop transmitting. Other remaining links in the network 
do not affect the time until the next state change since their 
interfering neighbours are transmitting. Moreover, since packet 



transmission times and back-off periods are exponentially 
distributed, the time until the next state change is the minimum 
of all of these associated exponential r.v.'s. In addition, since 
back-off timers and packet transmission times are independent 
of each other, given the state of system at time ?i£+X^=o _1 s i' 
the time until the next state change is the minimum of a 
collection of independent exponential r.v.'s. As a result, the 
r.v. Si is also an exponential r.v. with a rate that depends only 
on system history up to time ne+J^'jZo 1 s j- Since there are at 
most L links, for z > 1, the r.v. s$ becomes an exponential r.v. 
with maximum rate of Lz. Hence, independent of the system 
history up to time ne and the previous times where links states 
changed, r.v. s, is stochastically dominated by an exponential 
r.v. with rate Lz and we have that 

E Si {sj_i, Si_ 2 , si},H s (ne) > 
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Lz 



(288) 



where W s (t) is the history of the system up to and including 
time t. 

Define v(L) as 

z\f~L = o(zL), as L — > oo. 



v{L) 



(289) 



Multiplying both sides of ( |287| i by —v(L), taking exponential 
function of both sides, and then applying Markov inequality, 
we obtain 



cLr) 2 

P[^2si<e \% s {ne) 

i=l 



Efe-^ii 2 *. \H a (ne)] 



< 



(290) 



In the above and the following to simplify the notation, we 
have used u as a short notation for v(L). Using ( |288[ ), and 
noting that for an exponential r.v. s with mean E[s] 



E[e- Vs 

we can show that 



(1+vE[s])- 



E[e 



\U s {m)\ <{! + —) 

XjZ 



-cLt] 2 



(291) 



For a fixed z and sufficiently large L, by ( j289[ >, we have 
that 



(1 



-cLr}2 



= (i 



< 



v 

Lz' ' Lz 

Hence, by the above and ( |290[ ), for sufficiently large L. 

cLr) 2 



(292) 



P[^2s l <e \H(ne) 



i=l 



< e 



^(cr/2-zs) 



(293) 



Using ( |289| > and noting that by ( |45| > for a fixed 

0, 



lim ze 

L— >oo 



we can choose L sufficiently large so that 



cLr] 2 



P 



\ Si < e \H s {ne) 



»=i 



< e 



CVXT72 



(294) 



The discussion leading to ( |287[ ) states that in order to have 
(|285|) for some t € [ne, (n + l)e), the inequality in (|287[> 



should hold. Using the probability bound in the above, for 
sufficiently large L, we thus have that 



P 



3t e [ne, (n + l)e), s.t. 
\\9 L {t)-e L {ne)\\ > m \H s {ne) 



< e 



(295) 



Note that the above probability bound is independent of the 
system history up to and including time ne. Hence, considering 
the event that in none of the intervals [ne, (n + l)e), < 
n < re -1 — 1}, the inequality of ( |285| l holds for some t E 
[ne, (n + l)e), we have 



P 



sup sup \\0 L (t) - 9 L (ne)\\ < ?/ 2 

0<n<f-l ne<t<(n+l)e 



> [I - e 

Since e = [L~^-~^], it follows that 

(1-c - ) =1. 



lim 

L— >c 



(296) 



(297) 



Using this limit and the bound in ( |296| l, we obtain the 
statement in the lemma, as required. 

□ 

Lemma 13. Under the assumptions in Lemma [70] there exists 
a constant zq such that following holds for z > Zq. There exists 
a ti independent of z, where 

t <h< 2, 

such that x(i) uniquely exists over the interval [to,ti], and for 
all t G [to, t{\, we have x(i) G T> w . Moreover, we have 



< t£ i < xi(h) < c' xu2 < 0.5, 



(298) 



where c' Xl 1 and c' Xl 2 are constants independent of z, and for 

2 < i < i 



x l (h) = o(z- 1 ). 



In addition, 

x 2 {ti) 



3z 
2 

X3(tl) - — Cfl<5 x (*l) 

iz 

Xi{h) - - [xi(h) - c R 5 x (hf\ = 0(z-' A ). 



= 0(z-% 
= 0(z- 2 ). 



(299) 

(300) 
(301) 
(302) 



Proof. By the assumption in the lemma on x to = x(to), for 
large z, 



|x(to)|| < 1. 



(303) 



In the rest, since the lemma is stated for z > Zq, we assume z$ 
is chosen sufficiently large so that the above inequality holds 

for z > z . 

Consider the ODE of < | 1 1 0| > and suppose z > zq. Since the 
RHS of this ODE is a locally Lipschitz continuous function 
of x and continuously differentiable with respect to Xi's, 
1 < i < 4, we have that starting from the bounded x to , see 
P03|l, the solution exists and is unique until at least where we 



have ||x(f)|| > 2 for some t > t (see e.g., Theorem 1.9 and 
Theorem 1.13 in [45]). Define time tb to be the first time after 
t that 

||x(t 6 )|| > 2. 

By the definition of tb, x(t) uniquely exists over the interval 
[to,tb), and we have that 

||x(t)|| < 2, t G [t ,t b ). (304) 

If time tb does not exist, define tb = oo. By ( |303| ), we have 

> to. 



By the definition of /(x) in < | 1 1 0) > and the inequality ( |304| i, 
we also have that 



|/(x)| < cl 4x i, t G [t , 4), 



(305) 



for some constant c > 0. By the ODE of \\ 10[ ), we also have 
that for all t G [to, 4), 

x(t) = x(t )e A (*- t0 ) + / e A ( t -*')/(x(t'))^'- (306) 

J to 

Since all eigenvalues of matrix A are real, distinct, and non- 
positive, we can find a constant ca > independent of z such 
that for t' < t 



< Ca14x4- 



(307) 



Considering a time t G foi;i&), by ( |304| i-( (307l >, we have 



|x(t)| < 2c A l4xll4x4 + / 
J to 



f 


e A(t-t') 


>t 





|/(x(f))|df 



< 8ca14x1 + CCa(* — t )l4x4l4xl 

<4cc A (t-to + 2c- 1 )l 4 xi. 



(308) 



The above inequality states that each element of x(t) 
independent of z increases at most linearly with time t before 
reaching time %■ Therefore, by ( |303[ ), there must exist a time 
interval 

^0 = fob U] 



such that 

||x(t)|| <2, te/o 

where f j can be chosen independent of z such that 
to < U < 2. 

As a result, 

\ Xi {t)\ < 2, t G lo- 
Using the above inequality, we have that 

\Sl(t)\ <8, tG/ 

where 5 x (i) is defined in ( |109| l. 

By the ODE for x^if), as given by (jTTOjl, we have 

d 



dt 



X2 = Cr5^ — 3ZX2- 



(309) 
(310) 
(311) 
(312) 

(313) 
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Using the above ODE, ( |312[ ), and the assumption in the lemma 
that ( |215| l holds, i.e., that ^(to) = we can show that 



sup|a:2(t)| =0(z _1 ). 
tei 



(314) 



Similarly, we can show that 



sup \x 3 (t)\ = O^" 1 ), sup |x 4 (t)| = 0(z- 1 ). (315) 



Using pTT) , ( [3141 , and ([3T3J in the ODE for x x {t) given 
in d 1 1 0b, we have that 



sup 

tela 



dt 



Xi 



= 0(1), as z 



Hence, for any e\ > 0, we can find a sub-interval 

h = [*d,*<] da 

with 



(316) 



(317) 



to < t'i < U, 

and t\ independent of z such that 

\x x {t) - x x {to)\ < ei, t G Ji. 

In particular, since by assumption in the lemma 

< c Kli i < ii(i ) < c Xu 2 < 0.5, 

we can choose t\ sufficiently small such that for all t G 1\ 

0< < 4 lil <a:i(t)<<4 li2 <0.5 ) (318) 

where d x j and c^. 2 ^ constants independent of z. 

Using ( |318| l, the definition of <5 x (t) given in ( |109| l, the 
assumption in the lemma that x?,{to) = 0{z^ r ), and the ODE 
for X2(t), given in ( |3 1 3b , we can show that 



< x 2 (t), t e h, 
sup x 2 (t) = 0(z~ 1 ). 

Similarly, for 3 < i < 4, we obtain that 

< Xi(t), teh, 

sup Xi{t) = 0(z^ 1 ). 
teii 



(319) 



(320) 



Considering ( |318| >, ( |3 19b , and ( |320b , for z > z , we have 
found an interval I\ C Iq C [to,ifo) independent of z within 
which x(t) uniquely exists and is well-defined, i.e., we have 
that 



x(t)6D«, te/i. 



(321) 



We next show that ( |300| )-( [302] i hold. Here, we provide the 
proof for pOOl The proofs for pOl) and < |302} follow from 
similar lines. Consider the ODE for Xi(t) given in ( |3 1 3b , 
which can be used to express x 2 (t) as a function of cjiS^.(t) 
for t G ii as 

a; 2 (t) = x 2 (t )e- 3z * + /' e- 3 ^*-*')^^)^'. (32 2) 



to 



Recall that Jj = [to,^] where t\, t\ > t , is chosen 
independently of z. This allows us to choose a sufficiently 
small constant r\ > and choose zq sufficiently large so that 

for z > zq and k — y/z, we have 



t +V+- <t'i- 



z 

Consider a time t G (to + r) + | ,tj] C I\. Using ( |322] >, we 
have 

x 2 (t)= [ e- 3 ^-^c a 6l(t')dt' + e X2 (t), (324) 
where by P03[ ) and ( |312| i, we have 

e, 2 (t)-a; 2 (<o)e^ (t - to) 
< e -a»(t-to) + ^*( 



(323) 



e- 3z(t -* o) ). (325) 
By Mean- Value Theorem, we have that 



3z 



-3z(t-t') 



Cfl 4(t) + (t'-t) Cii -[4(t")] 



eft', 
(326) 



where t" is a function of i' and t, and we have that 

t" e(t',t) dhC lo- 
Using ( fT09) >, pl2) , and ( ffTB) , we have that 



sup 



4(0 



O(l), as z — > oo, 



(327) 



which can be used along with ( |326| > to show that 



sup 

te(to+i7+£>*'i] 

< 0(l)e 



3z 



-3A: 



■0(l)(z) 



(328) 



Since fc = y/z, combining ( |324| i, ( |325| >, and ( |328| >, we obtain 
that 



sup 

te(t +v+i,t'i\ 



0(z~ 



(329) 



Since t\ does not depend on z, we can choose 77 sufficiently 
small and zq sufficiently large such that for z > zq, we have 

k 



Hence, by ( |329| l, we have that ( 300| > holds for t\ as defined 
above. Similarly, we have that ( |301) > and ( |302[ > hold. Moreover, 
by the definition of t\, we have that for z > zq 

[t ,h] c h. 



Hence, by (f3T8]>-<|320|), we have (p98]l-(|299l hold. Moreover, 
since [io,ii] C I\, by ( |321) >, we have that for z > z , x(i) 
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is uniquely defined over the interval |io,iiJ, and that for all 
t G [tojii]> x (*) G D w , completing the proof. □ 

Lemma 14. Suppose y(t) and x^(t) are such that we have 
d 



dt 



y = Ay + /i(y) + e y , 



dt 



x h = Ax h + h(x h ), 



(330) 



with the initial condition at time t\ < r f/zaf 

Xfc(ti) = y(*i) = y tl 
/or some constant vector y tl , ant/ where 

sup ||e y (t)|| = 0(z _1 ). (331) 

te[ii,r] 

Moreover, suppose the function h(-) is a Lipschitz continu- 
ous function with the property that for any Xi and x 2 



ft-(xi) - /i(x 2 ) = L /l (xi,x 2 )[xi - x 2 ] 



(332) 



where L/j(x!,x 2 ) is a 4x4 matrix with elements 
[L/ l (xi, x 2 )]y such that 



(333) 
(334) 
(335) 

(336) 



where by ([333]), the matrix L/j(i) is such for the absolute value 
of its element [Lh(t)]ij, we have 



|[Lfe(xi,x 2 )]ij| < L h , 
for some constant Lh > 0. Let 

e h (t) = y(t)-x h (i). 
Wfe fnen have that 

sup ||e h (t)|| =0(^ _1 ). 
terti.ri 



Proof. Using ( |332| i and ( |334| i, we have that 

ft(y(t))-^(xh(t)) = L h (t)e h (*) 



|[M*)U < 



(337) 



Subtracting the ODE for x g (t) from that of y(t), as given 
in the lemma, we obtain 

d 



dt 



e h = Ae h + L h e h + e y = [A + L h ] e h + e y . (338) 



By the initial conditions in the lemma, we have 

Xft(ti) = y(*i), 

and hence, 



eh(h 



0. 



Using this initial condition for eh(t) and the ODE in < |338) >, 
we have 

e h (t)=V(t) f V(t') _1 eyW (339) 

where V(i) is the fundamental matrix solution for the ODE 
in ([338} and is given by 



V(t) = e #i [A+I * (r)l *'. 



(340) 



Plugging this into (|339|>, we have 



e h (t)= / eJ^ A+ ^ {r ^ dr e y (t')dt'. 



(341) 



Since all eigenvalues of matrix A are real, distinct, and non- 
positive, we can find a constant ca > independent of z such 
that for t' < t 



,(t-t')A 



< caXa 



(342) 



where the inequality is component-wise. Using the assumption 
in ( |33 1| >, we then have that 

f °ti Adr e y (t')dt' 



sup 

te[ti,Tl 



< sup 



/' 


e (*-*')A 


Itl 





\e y (t')\dt' 



< (r - h)0{z- l )U xl = 0(«- x )l 4 xi. (343) 
By ( |337| >, we also have that 

|Lh(t)| <L h U xi , (344) 
which can be used to show that independent of z 

SUp \ef*' Lh ( r ) <fr | < e It' L h 1 ix4dr _ e (T-t')L h Uxi 
*e[ti,r] 

<c /l e 4i "( T - t ')l 4x4 (345) 

where Ch > is a constant independent of z. Therefore, by 
( |331| ), we have that 



sup 

te[ti,r] 



<c h I < 



AL ^-^U^\e y {t>)\dt> 

-l^«h(r-t 1 ) 1 



= 0(z- ± )e^ T -^l ixl 
= 0(z- 1 )l 4xl . (346) 
Using P41| >, ( |343| >, and ( |346| >, we obtain that 

sup \e h (t)\ =0(z- 1 )l M , (347) 

*e[ti,r] 

which implies the statement in the lemma, completing the 
proof. 

□ 

Lemma 15. Suppose for some t\ < t, we have y(ii) G T> w 
such that 

0<c; iiX <yi(h) <4x,2 < 0.5, (348) 

where c yi 1 and c' yi 2 are constants independent of z. In 
addition, suppose for I = {t{\, we have 

(349) 
(350) 
(351) 
(352) 



sup max [y 2 (t),y 3 (*),2/4 (*)] =0(2 1 ), 



sup 

te/ 



sup 

te/ 

sup 

te/ 



Ifc(t) - 2jfc(t) =0(0- 
y4(t)--[tfi(*)-c fl iy(t) 3 ] = 0(z~ 2 ) 



We f/ien have that y(t) G T> w , for t G [ti, r], ant/ 
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sup 

t£[ti,r] 



df 



0(1)1, 



(353) 



Moreover, we have that ( [349] > - <|352]> hold for 

I=[h,T). 

Proof. Using the ODE for yi(t) in ( |1 13[ ), we obtain 

Vl (t) = 0.5 - a' (I + /3' (t-h))-^ 

where 

a' = 0.5-y 1 (t 1 ), 

and 



/3' 



3 v ; 



(354) 



(355) 



(356) 



Having determined yi(t), we then can use ( |1 1 3| > and yi(t) 
to find y 2 (i), j/3 (t), and yi(t). Since by ( |354| >, 



< j/i(t) < 0.5, 



(357) 



we can use the constraints on the initial conditions given by 
( |349| l at time i = t% to show that j/2(f), 2/3 (i), and j/4 (t) are 
all positive for t G [ti,r], and that ( |349] > holds for / = [ti, r]. 
Therefore, y(i) G for t G [ti,r]. Having (357} and that 
((349) holds for I = [ti,r], it follows from < fTT3> that 



sup 

te\t 1: T] 



df 



= 0(l)l 4 xl. 



(358) 



We next prove that ( |352[ ) holds for / = [ti,r]. The proof 
for ( |350[ ) and ( |351| l follows from similar steps. By ( |358| l, we 
have 



sup 

te[*i,T] 



df 



yi(t) = O(l), as z —> 00. 



(359) 



Now consider the ODE for 2/4 (i) given by ( |1 13[ ). We have that 

j t Vi = [2/1 - Cfi[0.5 - t/i] 3 ] - zy A . 
Define yd(t) as 

y<j(*) = 2/1 - ch[o.5 - z/i] 3 . 

By ((357) and (|359|, we have that 
d 



(360) 



(361) 



sup 

*e[ti,r] 



(7f 



O(l), as z — >• 00. 



(362) 



(363) 



By Mean- Value Theorem, for f G [ti, r], we have 

Vd(t) = + (* - *i)4;i/d(*i) 

where ti is a function of t\ and i such that ti G (ti, i). 
Knowing 2/1 (£) and using P60[ ), we can find 2/4 (i) as 

y 4 (t) = y4(ti)e-* ( *- t0 + / e-^"*'^ (364) 



To obtain the statement in the lemma, we use ( |364| i. By ( |363[ ) 
and ( |364| i, we have 

y A {t) = y^e-^-^ + f e - z ^y d (t)dt' + e ViA (t) 

Jti 

(365) 



where 



e s/4,i(*) = / ' 
Iti 



s(t_t,) (*' -t)^y d (t")dt', (366) 
at 



in which t" is a function of t' and t such that t" € (t',t). 
Using (|362), we have that 



sup |e y4! i(t)| = (9(z 2 ). 
te[ti,T] 



(367) 



We also can write 



-z(t-tx) 



-z(t-t') 



Vd(t)dt' 



ti 



= e -*(*-*0 [j, 4 (t x ) - i w (t)] + i w (t). (368) 
L z z 



Moreover, using ( [362} and ( 563] ), and noting that xe xz < 
z _1 e , we have that 



sup 

t€[tl,T] 



,-*(*-*!) I W ( t )_ e -*(*-*0 I 



yd(ii) 



< sup e" z(t - tl) -(i-ti)0(l) = 0(z~ 2 ). 
te[ti,r] z 

(369) 



Using ( [365) , p67] ), ((368), and ggg) , we have 



sup 

te[ti,r] 



y4(i) Z/d(*) 

z 



< sup e 

te[tl,r] 



-a(t-ti) 



J/4(*i) 2/d(*i) 



0(z" 



(370) 



Using the above and the assumption in the lemma given in 
(j352|) for J = {ti}, we then have that 



sup 

tefti.rl 



2/4 0) - 



as required. 



0(z~ 2 ), (371) 



□ 



Lemma 16. Consider the lattice Gl, attempt rate z > 0, one/ 
a g/ve« f;me i, < t < oo. We /;ave, w.p.l, 

£cec[" d) (M) l ^ 
where T)(L) > and 



>S L (t)~r,(L)~59 Wl (t) 



lim ?7(X) = 0, 

I— foo 

one/ 0L,h(t) is the density of links that are inactive and sense 
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Lemma 17. Consider the lattice interference graph Gl, 
attempt rate z > 0, anJ a time t where < t < oo. Recall 
that 0L,{t) represent both density and the set of active links at 
time t. Let 

be the total area of the lattice not covered by the union of the 
coverage areas of active links at time t, where A\ is defined 
in Appendix^ Then, w.p.l, 

A c L (t) 



2L 



-s L (t) 



< V (L) 



(372) 



where rj(L) only depends on L, and 
lim T](L) = 0. 

Proof. By the definition of A L (t), the area covered by active 
links at time t is n 2 — A c L (t). Separating the contribution of 
active links on the boundary 8Gl of the lattice Gl from those 
of active links inside the lattice with Ai = 2, we obtain 



A c L (t) 



ie(dG L ne L (t)) k l^G L nlee L (t) 



0(VL)+2{1: I £dG L r\l e0 L (t)}\ 



(373) 



The term 0{yL) accounts for the contribution of links in 
8Gl, whose total number is less than AyL, and whose 
coverage area is less than 2. Similarly, we have that 

e L {t) = ^\{H B(G L ) n i e e L (t)} | + 0(1/ VI) (374) 

Dividing ( |373| ) by 2L, and using fl8] l and ( |374| i, we have 
1 



2L 



(n 2 -Al(t))-6 L (t) 

±-(n 2 Al(t)) - 0.5 + 6 L (t)\ = 0(1/VL). 

(375) 

Noting that n 2 /L = 1 - 0(1/ VZ) since L = (n + l) 2 , from 
the above, we obtain 

A c L (t) 



2L 



8 L (t) =0(1/VL), 



which implies the statement in the lemma, as required. □ 

Lemma 18. Consider the lattice Gl, attempt rate z > 0, and 
a given time t, < t < oo. We have, w.p.l, 

A%(t) <V2 J2 £ ( C ) + 10X ^(*) 



ceei" d) (M) 



where A c L (t) is defined in Lemma 17 C L nd \t, z) is the set of 
non-dominating clusters at time t as defined in Section \VIII\ 
and 0L,h(t) is the density of links that are inactive and sense 



the channel as idle at time t, as defined in Appendix II-B.2 the channel as idle at time t, as defined in Appendix \II-B.2 



Proof. The lemma follows from Lemma 17 and Lemma 18 



□ 



Proof. Define the set £^(t) = {l hfi , 1 <i< i h } to be the 
set of all inactive links that sense the channel as idle at time 
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t. By definition, LOl h(t) is the total number of these links, 
and hence 

i h = L6 L<h (t). 

Consider the following process. First, suppose ih > 0; we 
later also consider the case of ih — 0. Start from the first link 
Ih i in L^ h \t) and make it active, which results to a new set 
of active links. We then consider the next link lh,2- If this link 
does not cause interference to the new set of active links, we 
make link 1^2 active. We continue the same process for all 
i = 3, ih- By this process, 1) we reduce the uncovered area 
A c L (t) by at most 

2L6 L , h (t) (376) 

since link coverage area is at most 2, 2) we might have a new 
set of clusters, which we denote by 

C' L (t,z) = {C' i ,l<i<i' max h 

and 3) the length of cluster boundaries changes at most by 
4\/2L9l h(t) since the coverage area of each newly added 
link can change or contribute at most 4\/2 to the length of 
cluster boundaries. Therefore, 



t(C')- £ ( C )| <4V2Le L , h (t) (377) 

cec[ s) (M) cec { l ld) (t,z) 

where C ( L s) (t, z) is a subset of clusters in C' L (t, z) that have 
the same type as the non-dominating clusters in C^ d \t,z). 
For the case where ih = 0, we have that 6L.h(t) = 0, and 
the above inequality trivially holds by letting the new set of 
clusters to be the same as the original set of clusters at time 
t, i.e, by letting C' L {t, z) = C L (t, z). 

After the above process, we cannot add any further links 
to the set of active links, and the remaining uncovered area 
will be the union of areas around new cluster boundaries. For 
instance, consider the uncovered area between the two clusters 
in Fig. [9] By inspection, we find that this area is less than 

/2 

c u \m (378) 

where c u is a constant such that < c u < 2, and C can be 
either of the clusters in the figure. 

Considering the area removed in the defined process and 
the area between new cluster boundaries along with ( |376| > and 
( |378| ), we have that 

A c L (t)<V2 Y e{C') + 2L9 L . h (t). 

cec^\t, z ) 



Using this inequality and ( |377[ ), we obtain the inequality in 
lemma. □ 



